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ABSTRACT 

We present broadband (radio, optical, and X-ray) light curves and spectra of the afterglows of 
four long-duration gamma-ray bursts (GRBs 090323, 090328, 090902B, and 090926A) detected by 
the Gamma-Ray Burst Monitor (GBM) and Large Area Telescope (LAT) instruments on the Fermi 
satellite. With its wide spectral bandpass, extending to GeV energies, Fermi is sensitive to GRBs with 
very large isotropic energy releases (lO^** erg). Although rare, these events are particularly important 
for testing GRB central-engine models. When combined with spectroscopic redshifts, our afterglow 
data for these four events are able to constrain jet collimation angles, the density structure of the 
circumburst medium, and both the true radiated energy release and the kinetic energy of the outflows. 
In agreement with our earlier work, we find that the relativistic energy budget of at least one of these 
events (GRB 090926A) exceeds the canonical value of 10^^ erg by an order of magnitude. Such energies 
pose a severe challenge for models in which the GRB is powered by a magnetar or neutrino-driven 
collapsar, but remain compatible with theoretical expectations for magneto-hydrodynamical coUapsar 
models (e.g., the Blandford-Znajek mechanism). Our jet opening angles (9) are similar to those found 
for pre- Fermi GRBs, but the large initial Lorentz factors (Fq) inferred from the detection of GeV 
photons imply ^Fq « 70-90, values which are above those predicted in magnetohydrodynamic models 
of jet acceleration. Finally, we find that these Fermi-h AT events preferentially occur in a low-density 
circumburst environment, and we speculate that this might result from the lower mass- loss rates of 
their lower-metallicity progenitor stars. Future studies of Fermi-hAT afterglows in the radio with 
the order-of-magnitude improvement in sensitivity offered by the Extended Very Large Array should 
definitively establish the relativistic energy budgets of these events. 

Subject headings: cosmology: observations - gamma rays: bursts - radio continuum: general 
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Long-duration gamma-ray bursts (GRBfE!l), like 
hydrogen-deficient Type Ib/c supernovae (SNe Ib/c), re- 
sult from the gravitational collapse of the evolved core of 
a massive star. The main characteristic that sets GRBs 
apart from other SNe is that a substantial fraction of the 
energy of the explosion is coupled to relativistic ejecta. 
A compact central engine is responsible for accelerating 
and collimating the se jet-like outfiows and driving the SN 
explosions ( WooslcY_fcBlooin|[20p IGehrels et al][200l 
ISoderberFet al.ll2010D . The precise nature of the cen- 
tral engine which powers GRB-SNe, however, remains 
an open question. 

Motivated by empirical constraints, all viable central- 
engine models for long-duration GRBs share some com- 
mon characteristics (e.g.. [Piranll2005t) . They must pro- 
duce a collimated outflow with an initial Lorentz factor 
(Fq) of a few hundred on observed time scales of 10-100 s, 
with luminosities and kinetic energies of order lO^'' erg 
s~^ and 10^^ erg, respectively. Leading models include 
the "collapsar" model in which a relativistic jet is pro- 
duced from a r o tating black hole/accretion d isk system 
(jWooslevI [l99l iMacFadven fc WooslevI 11999'). and the 
"magnetar" model in which the rapid energy loss from a 
newly born millisecond neutron star (formed either from 

Throughout this work, we shall refer exclusively to long- 
duration GRBs (i.e., those apparently with massive-star progen- 
itors) unless explicitly stated otherwise. 
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the gravitational collapse of a massive star or from ac- 
creting or coalescing white dwarfs) with a 10^^ G mag- 
netic fiel d drives a Poynting flux-dominated relativistic 
outflow (|Usovlll992[ ). 

These and other more exotic models for GRB cen- 
tral engines are highly constrained by their energet- 
ics. The prompt high-energy emission, when com- 
bined with a spectroscopically determined redshift (and 
hence distance measurement), yields the isotropic ra- 
diated gamma-ray energy (Ej^iso)- Well-sampled af- 
terglow observations allow both a measurement of the 
degree of coUimation (and hence the true beaming- 
corrected energy release in the prompt emission, E^j 
and the kinetic energy remaining in the shock that 
powers the broadband afterglow emission (-Eke)- Such 
measurements, made nearly a decade ago, pointed to 
a total relativistic energy yield (ij^rei ~ -\- Eke) 
of ~ 10^^ erg dlVail et alJ 120011: iPanaitescu fc Kumail 
200Ti iFreedman fc WaxmanI l2001t iBloom et alj 120031: 
Berger et al.ll2003aD . 

Since that time, there has been growing evidence for 
a considerable range in the relativistic energy scale i?roi, 
suggesting either a diversity in central engines or their 
properties. Most notably, a population of nearby (red- 
shift z < 0.1) subenergetic long-durati on GRBs have 
been identifled (jBloom et al. ; 2003; Sodcrb erg et al.ll2004 
120061) . They too are associated with SNe Ib/c, but their 
relativistic energy release is a factor of 100 less than that 
of typical cosmological GRBs and their outflows are sig- 
nificantly less coUimated (quasi-spherical). Since they 
can only be detected at low redshifts where the com- 
parative volume for discovery is low, they are small in 
total number. But their volumetric rate is inferred to be 
10-100 times la rger than that of the more distant long- 
duration GRBs (fSoderberg et al.l[200llCobb et alll2006l : 
iLiang et al.ll2007[ r 

More recently, evidence has been growing for a class of 
GRBs whose total relativistic energy release is at least an 
order of magnitude above the canonical value of 10^^ erg 
(e.g., .Cenko et al. 2010. and references therein). Unlike 
subluminous events, the total energy budget of these 
hyper-energetic events poses a significant challenge for 
some progenitor models. In particular, models in which 
the GRB is powered by a magnetar or a neutrino-driven 
collapsar are strongly disfavored. On the other hand, 
collapsars driven by magneto-hydrodynamical (MHD) 
processes, such as the B landford-Znajek mechanism 
(jBlandford fc Znaield Il977| ) , can naturally accomodate 
energy budgets as large as 10^'^ erg. 

Unfortunately, it has been rather difficult to con- 
strain the beaming-corrected energetics fo r the hundreds 
of GR Bs detected by the Swift satellite (jGehrels et al.l 
[200l . The reasons for this difficulty are now 
largely understood. First, the relatively narrow en- 
ergy bandpass (15-150 keV) can miss entirely the 
peak of the gamma-ray spectrum, making estimates 
of E^^iso highly uncertain. Second, there has been 
a dearth of mea s urements of jet opening angles (e.g., 
Panaitescul 12 0071: iKocevski fc Butle'il I2Q08I: ILiang et alJ 
2001 " Racusi n et al. I I2009D and well-sampled multi- 
wavelength GRB afterglows (used to derived the after- 
glow kinetic energy Eke)- Swift GRBs are on av erage 
more than twice as distant (| Jakobsson et al.l 120061 ) and 




Redshift 

Fig. 1. — Prompt isotropic gamma-ray energy release (-E-^jso) of 
GRBs. With its soft, narrow bandpass (15-150 keV), Swift typi- 
cally selects events with smaller isotropic energy release but larger 
opening angles than previous missions, which t riggered predomi- 
nantly in the MeV bandpass l|Perna et al.l[20i3 ). GRBs detected 
at GeV energies with the Fermi-LAT all fall at the brightest end 
of the isotropic energy distribution, and must therefore be highly 
coUimated to achieve a canonical beaming-corrected energy re- 
leaae of ~ tO^^ srg . References: pre-5'u'i/f: Amati ( 2006); Swift: 
IButler et all 1 120071) : Fermi-LAT: IGreiner "eFaLl lf2009l ), this work. 

therefore significantly fainter (^ 1.5 mag in the optical; 
iBerger et alll2005t iKann et al.ll2007l) than GRBs in pre- 
vious samples. In large part this is due to selection ef- 
fects: a combination of bandpass and sensitivity from 
Swift has preferentially selected the faint end of the lu- 
minosity function — GRBs wit h low isotropic ene rgy re- 
lease but large opening angles (jPerna et al.ll2003t) . 

With its nearly seven decades in energy coverage 
(lOkeV - 100 GeV), Fermi can provide unparalleled con- 
straints on this subsample of the most luminous events. 
In light of the empirical relation between the peak 
energy of the gamma-ray spectrum and the isotropic 
gamma-r ay energy rel ease (the Ep-Ej^iso, or "Amati" 
relation; [Am atil 120061 ) . MeV/GeV events detected by 
either t he Gamma-Ray Bu rst Monitor (GBM, 8keV - 
40 MeV; iMeegan et al. 20091 or the Large Area Telescope 
(LAT, 20 MeV - 300 GeV; Atwood et al l [2009, ) onboard 
Fermi preferentially select a sample of GRBs with large 
isotropic energy release (Figure [TJ. IIigh-£^.^4so events 
also h ave brighter X-ray and optical afterglows on aver- 
age ( Nvsewander et aTTl2009f l. Follow-up afterglow ob- 
servations can then determine whether these GRBs are 
highly beamed events {9 < 2°) with a typical energy re- 
lease or true hyper-energetic GRBs. 

The Fermi-hAT offers a further advantage over previ- 
ous GRB missions sensitive only at MeV and keV ener- 
gies by providing strict constraints on the initial Lorentz 
factor of the relativistic outfiow. To avoid e"*" — e~ 
pair production (and the accompanying thermal spec- 
trum), the GRB jet must be moving towards the ob- 
server wi th ultra-relativistic sp eeds (the "compactness" 
problem; iCavallo fc Reeslll978D . The higher the energy 
of the most energetic photon detected from a GRB, the 
more strict the lower limit on the outflow Lorentz fac- 
tor will be. Combining the Lorentz factor limits for the 
most relativistic GRBs with inferred jet opening angles 
from broadband afterglow models can provide critical di- 
agnostics of the jet acceleration mechanism. 

Here we report on broadband (radio, optical, and 
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Fig. 2. — The broadband radio (blue), optical (red), and X-ray (black) light curve of GRB 090323. The best-fit model is plotted in solid 
lines (see Table |4] for parameters). The identical model parameters for an isotropic explosion are plotted as the dashed lines. The strength 
of the possible modulation of the radio afterglow caused by interstellar scintillation (e.g.. Frail ct al. 2000a) is indicated by the light-blue 
shaded region. The model provides a reasonable fit in all bandpasses. It is clear that any jet break must occur at t > 10 days, although 
the upper bound on the jet break time is only weakly constrained. 



X-ray) observations of four long-duration GRBs de- 
tected by the Fermi-LAT at GeV energies: GRBs 090323, 
090328, 090902B, and 090926A. For each event we con- 
struct afterglow models to constrain the collimation 
and beaming-corrected energetics, and we compare these 
LAT events with previous GRBs detected at other en- 
ergies (i.e., keV energies from Swift, and MeV energies 
from pre- Swift satellites). For three of these GRBs, we 
also present the optical spectra used to determine the af- 
terglow redshift. A more thorough analysis of the host- 
galaxy properties of these events will be presented in a 
forthcoming work. 

Throughout this paper, we adopt a standard ACDM 
cosmology with Hq — 7 1 km s~^ Mpc""^, fl^ — 0.27, and 
f^A = 1 - Qni = 0.73 (ISoergel et all [2001 . We define 
the flux-density power-law temporal a nd spectral deca y 
indices a and /3 as /j, oc t~°'iy~^ (e-g-, iSari et al]|1998[ ). 
All quoted uncertainties are la (68%) confidence inter- 
vals unless otherwise noted. 

2. OBSERVATIONS 

2.1. GRB 090323 

2.1.1. High-Energy Properties 

GRB 090323 was detected b y the Fermi GBM at 
00:02:42.63 on 23 March 2009 ()Ohno et al.l [200l UT 
dates are used throughout this work). In the 8keV to 
40 MeV bandpass of the GBM, the light curve was mul- 
tipeaked with a duratio£3 of tgo w 150 s. GRB 090323 
was also detected at MeV energies by several of the 
satellites comprisin g the Inter-Planetary Network (IPN; 
iHurlev et al.lll999() , including the Konus instrument on 
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the W ind satellite (jHurlev et al.ll2009t iGolenetskii et al.l 
[2009^ . 

Examining only t he fi rst 70s of GBM dat^, 
Ivan der Horst &: XinI (|2009[ ) find that the prompt spec- 
trum is well described by a power law having an expo- 
nential cutoff at high energies with a — —0.89 ±0.03 and 
Ep = 697±51 keV. The resulting fiuence in the 8-10^ keV 
(observer frame) bandpass is fj = (l.OOiO.Ol) x 10""* erg 
cm~^. The Konus- Wmc? instrument was able to mea- 
sure the high-energy spectrum over the ent ire duration 
of the MeV emission. Fitt ing a Band model ( Band et aD 
119931 ) to their spectrum, IGolenetskii et al.l (12009d) find 
a reasonable fit with a = — 0.96lo.o9j 
Ep = 416t^^keV, and = 2.02;j;;^^ x IQ-'^ 
(20-10^ keV observer-frame bandpass). 

In order to facilitate direct comparisons between 
events, we must transform this high-energy fiuence to 
a common (rest-fr ame) bandpass (i.e., a "k" -correction; 
iBloom et al.l I200H ) . Here we adopt the rest-frame 1- 
lO^keV bandpass, as this encompass es the full range o f 
peak energies observed from GRBs (I Band et al.l Il993| ) 
without (for most pre-Fermi satellites) requiring large 
extrapolations outside the observed bandpass. At z = 
3.568 (§ I2.1.3p . The Konus- K^md measurement corre- 
sponds to a prompt fiuence of fj = (1.50 ± 0.20) x 
10~'*ergcm^^ in the 1-10^ keV rest-frame bandpass. 
Given the limited temporal coverage of the GBM spec- 
trum, wc shall adopt this value for the prompt fiuence of 
GRB 090323 for the remainder of this work. 

In addition to the detection at MeV energies by the 
GBM, GRB 090323 was also detected at GeV energies 



^® It is customary to report GRB durations measured as the time 
between the arrival of 5% and 95% of the background-subtracted 
fiuence. This quantity is referred to as tgo- 



^"^ After this point, the Fermi spacecraft slewed to position the 
GRB at the center of the LAT field of view, causing rapid changes 
in the GBM background level. 
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Fig. 3. — GMOS-S optical spectrum of the afterglow of GRB 090323. The broad absorption feature at A 5570 A results from Lyo 
in the GRB host galaxy. We identify a number of strong, narrow absorption features redward of hya from the GRB host galaxy at a 
common redshift of 3.568 it 0.004 (black annotations). The emission blueward of hya is strongly affected by the Lya forest. We identify 
two additional intervening absorbers, based on Mg II AA2796, 2803 (z = 2.101; blue annotations) and C IV AA1548, 1550 {z = 3.379; red 
annotations). Areas of the spectrum affected by telluric absorption in the atmosphere of the Earth are marked with crossed circles. Because 
of the increased noise in the NIR due to fringing, we have cut off the plot at A = 9000 A. 

filters, implied a spectral steepening around the observed 
g' band. Associating this br eak with abso rption from 
Lya in the GRB host galaxy, lUpdike et all (r2009b) de- 
rive a photometric redshift of 4.0 ± 0.3 for GRB 090323. 

We began observations of the afterglow of GRB 090323 
with the automated Pal omar 60 inch (1.5 m) telescope 
(P60; 'Cenko et al.' 2006") beginning on 2009 March 24 
(Cenko & Perley 2009). Images were obtained in the 



by the Fermi LAT. Li ke several previous GRBs observed 
at Ge V energies (e.g., Hurlev et al. 1994: Gon zalez et al.l 
l2?)03l [Giuhani et al. 2008: Abdo et al. 2009bi ). the GeV 
emission began several seconds after the MeV emission 
is detected, and remains significant long 1000 s) af- 
ter the MeV comp onent has faded beyond detectability 
(jOhno et al.ll2009f ). The highest energy photon from the 
direction of GRB 90323 had an energy of £; « 7.5 GeV 
(jPiron et al.ll2009l) . However, this photon arrived after 
the prompt MeV emission ended, and so its origin (i.e., 
prompt emission or afterglow, or, in physical terms, in- 
ternal or external shock-related) is somewhat uncertain. 
The highest energy photon detected by the LAT dur- 
ing the formal MeV igo measure ment had an energy of 
E w 500 keV (at t « 90s: IPiron et al. 2009). 

2.1.2. Afterglow Observations 

The Swift satellite began observations of the field of 
GRB 090323 with the onboard X-ray Telescope (XRT; 
iBurrows et all l2005b[) at 19 :27 on 23 March 200 9 (- 
19.4hours after the burst; iKennea et al.l l2009aD . A 
candidate afterglow was promptly identified at a = 
12h42'"50.''26, S = +17°03'14'.'2, with a 90% contain- 
ment radius of 2'.'7 (J2000.0). The XRT continued to 
monitor the evolution of the fading counterpart for the 
following two weeks. We have obtained the XRT light 
curve from the online compilation of N.R.BH the re- 
sulting evolution is plotted in Figure [2] 

The optical afterglow of GRB 090323 was discovere d 
shortly thereafter by GROND (Updike et al.l I2009b[ ). 
The optical/NIR spectral energy distribution, con- 
structed from simultaneous imaging in the g'r'i'z' JHK 

http: //astro . berkeley . edu/^nat/swif t; see 

IButler fc Kocevskil |[200W) for details. 



Sloan r' and i' filters and individual frames were au- 
tomatically reduced using our custom IRArF^ software 
pipeline. To increase the signal-to- noise ratio (S/N), in- 
dividual frames were astrometrically aligned using the 
Scamp software package and coadded using Swarjo- We 
used aperture photometry to extract the flux of the af- 
terglow from these coadded frames with the aperture ra- 
dius roughly matched to the full width at half-maximum 
intensity (FWHM) of the point-spread function (PSF). 
Aperture magnitudes were then calibrated relative to 
field source s from the S loan Digital Sky Survey Data 
Release 7 ( Abazaiian et a l. 2009) . Imaging continued on 
subsequent nights with P60 until the afterglow was below 
our detection threshold. The results of this monitoring 
campaign, uncorrected for for eground Galactic ex tinc- 
tion {E[B -V]= 0.025 mag; ISchlegel et a"LllT998l ). are 
presented in Table [HI 

We obtained additional imaging of the field of 
GRB 090323 with the two Gemin i Multi-Obje ct Spectro- 
graphs (GMOS-N and GMOS-S: iHook et al.l l2004'). Im- 
ages were taken at both Gemini North and Gemini South 
in the Sloan r' and i' filters, and they were reduced using 

IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association for Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

See http://astromatic.iap.fr. 



Afterglow Observations of Fermi-LAT GRBs 



5 



TABLE 1 
Radio Observations of GRB 090323 



Date 


At 


V 




Facility 


fUT) 


V J / 


(GHz) 


Vr^" J ) 




Mar 26.38 


3.38 


8.46 


27 ± 38 


VLA 


Mar 27.38 


4.38 


8.46 


225 ± 35 


VLA 


Mar 27.99 


4.99 


4.9 


105 ± 24 


WSRT'' 


Mar 28.43 


5.43 


8.46 


100 ± 40 


VLA 


Mar 29.16 


6.16 


8.46 


157 ± 31 


VLA 


Mar 30.18 


7.18 


8.46 


219 ± 39 


VLA 


Mar 31.32 


8.32 


8.46 


281 ± 38 


VLA 


Apr 1.30 


9.30 


8.46 


164 ± 35 


VLA 


Apr 3.29 


11.29 


8.46 


166 ± 27 


VLA 


Apr 3.28 


11.28 


4.86 


110 ± 45 


VLA 


Apr 4.41 


12.41 


8.46 


183 ± 35 


VLA 


Apr 5.14 


13.14 


8.46 


123 ± 29 


VLA 


Apr 6.28 


14.28 


8.46 


312 ± 27 


VLA 


Apr 7.42 


15.42 


8.46 


43 ± 27 


VLA 


Apr 9.44 


17.44 


8.46 


127 ± 30 


VLA 


Apr 10.42 


18.42 


8.46 


295 ± 27 


VLA 


Apr 14.09 


22.09 


8.46 


178 ± 56 


VLA 


Apr 17.04 


25.04 


8.46 


167 ± 33 


VLA 


Apr 25.46 


33.46 


8.46 


78 ± 32 


VLA 


May 3.41 


41.41 


8.46 


77 ± 31 


VLA 


May 11.35 


49.35 


8.46 


61 ± 29 


VLA 


Aug 1.18 


131.18 


8.46 


-13 ±27 


VLA 



Reference: Ivan der HorsH I I2009I) . 



the IRAF gemini package. Photometry was performed 
with the same methodology used for the P60 imaging of 
the field, and the resulting measurements are shown in 
Tabled 

Additionally, we have compiled optical and NIR mea- 
surements of the afterglow of GRB 090323 from the 
GCINEB Circulars and included these in Table [S The 
majority of the late-time measurements (of most inter- 
est to us for modeling purposes) were obtained in the R 
band. We note that most of these measurements were 
calibrated with respect to a sin gle object ( 1070-0238439) 
froi n the USNO -Bl catalog (M onet et aLll2003| ). follow- 
ing |Kann|er^ (|20d 9tf). Using the SDSS o bservations 
of this field and the filter transformations of iJordi et al.l 
|2P06), we find an i?-band magnitude for this source 
of i? = 17.31 ma g (compared to an assumed value of 
17.15 mag from K ann et ani2009bl ). We therefore offset 
all reported i?-band magnitudes from the GCN Circulars 
by -f0.16mag. 

Finally, we observed the afterglow of GRB 090323 with 
the Very Large Arraj0 (VLA) beginning a few days af- 
ter the initial burst trigger. The flux-density scale was 
tied to 3C 286 or 3C 147 and the phase was measured by 
switching between the GRB and a nearby bright point- 
source calibrator. To maximize sensitivity, the full VLA 
continuum bandwidth (100 MHz) was recorded in two 
50 MHz bands. Data reduction was carried out following 
standard practice in the AIPS software package. 

Our initial VLA detec tion of GRB 090323 was reported 
bv lHarrison et al.l (|2009i). GRB 090323 was also detected 
at radio wavelengt hs by the Westerbork Synthesis Radio 
Telescope (WSRT; Ivan der HorsH 120091 ) . The full set of 
VLA measurements is listed in Table [T] In order to im- 
prove the S/N and to reduce the modulation of the light 

■^^ See http : //gen . gsf c . nasa . gov/gcn3_archlve . html . 

The Very Large Array is operated by the National Radio As- 
tronomy Observatory, a facihty of the National Science Foundation 
operated under cooperative agreement by Associated Universities, 
Inc. 



curve caused by interstellar scintillation (e.g.. lFrail et al.l 
l2000al) . we binned the data from adjacent epochs. These 
binned points were used for our afterglow modeling (§[3]) 
and are plotted in Figure [2] 

2.1.3. Optical Spectroscopy 

We began spectroscopic observations of the optical af- 
terglow of GRB 090323 with GMOS-S o n 24 March 2009 
(^ 2 9.93 hours after the GBM trigger; iChornock et alj 
|2009|) . We first obtained 2 x 600 s spectra with the B600 
grating and a central wavelength of 6000 A, providing 
coverage over the range ~ 4500-7500 A with a resolu- 
tion of 2.5 A. Immediately following these exposures, we 
obtained 2 x 600 s spectra with the R400 grating and 
a central wavelength of 8000 A, providing coverage of 
~ 6000-10000 A with a resolution of 5.0 A. 

The spectra were reduced in a standard manner using 
routines from the IRAF ge mini and specred packages 
(see, e.g.. lCenko et al.ll2008l and references therein for de- 
tails). Wavelength calibration was performed relative to 
CuAr lamps and then adjusted based on measured night- 
sky emission lines. The resulting root-mean square wave- 
length uncertainty was < 0.25 A for the spectra with the 
B600 grating and < 0.40 A for the spectra taken with the 
R400 grating. Telluric features were removed using the 
continuum from well-exposed spectrophotometric stan- 
dard stars (e.g., iWade fc Horng 119881 : iMatheson et al.l 
[2000). 

Flux calibration was performed relative to 
the spectrophoto m etric standard star LTT 7379 
(|Stone fc BaldwinI [1985 Bald win fc Stone 1984). 
We caution, however, that the standard-star observa- 
tions were conducted on different nights from the GRB 
observations (31 August 2009 for the B600 grating; 4 
August 2009 for the R400 grating), so the absolute 
flux calibration is somewhat uncertain (estimated to be 
- 30%). 

The resulting spectrum of the afterglow of GRB 090323 
is shown in Figure |3l The broad absorption feature at 
A w 5570 A is produced by Lya at z w 3.6. The spectrum 
blueward of this wavelength is dominated by absorption 
from the Lya forest, while the lack of Lya forest features 
redward of this transition indicates that it corresponds 
to the redshift of the GRB host galaxy. Unfortunately, 
the CMOS CCD chip gap over the range 5500-5525 A 
precludes an accurate measurement of the Lya profile 
and hence a determination of the H I column density. 

Redward of Lya, we find a series of strong (rest-frame 
equivalent width Wq > 1 A) absorption features super- 
posed on a relatively flat, featureless continuum. In par- 
ticular, we identify Si II A1260, Si 11* A1264, O I A1302, 
Si II A1304, Si 11* A1309, C II A1334, C II* A1335, Si IV 
AA1393, 1402, Si II A1526, Si II* A1533, and C IV AA1548, 
1550 at 2 = 3.568 ± 0.004. Excited fine structure lines 
are indicative of UV pumping and further suggest these 
features arise in the host galaxy of GRB 090323. The 
strongest features exhibit a complex velocity structure 
that is only marginally resolved at the resolution of our 
B600 spectra. We also detect intervening absorption sys- 
tems of C IV AA1548, 1550 at z = 3.379 and Mg II 
AA2796, 2803 at z = 2.101. A more thorough analy- 
sis of the host-galaxy properties of GRB 090323 will be 
presented in a subsequent work. 



6 



Cenko et al. 




10° 10' 10= 

Time Since Burst (d) 



Fig. 4.— The broadband radio (blue), UV (red), and X-ray (black) light curve of GRB 090328. The best-fit model is plotted in solid 
lines (see Table [5] for parameters). The identical model parameters for an isotropic explosion are plotted as the dashed lines. The radio 
light curve is not very well fit at early times, although it likely suffers from strong interstellar scintillation (light-blue shaded region). In 
this case a jet break is required by the data to fall at t > fOdays. 



At z — 3.568, the isotropic prompt energy release 
from GRB 090323 in the rest-frame 1-10'' keV bandpass 
is E^Aso = (3.99±0 53) x 10^-^ erg. Using the formulation 
of iLithwick fc Saril (|2001[ ) , the lower limit on the outflow 
Lorentz factor (assuming a nonthermal spectrum up to 
Eohs « 500 keV) is Tq > 600. 

2.2. GRB 090328 
2.2.1. High-Energy Properties 

GRB 090328 tri ggered the Fermi GB M at 09:36:46 on 
28 March 2009 (,McEnerv et all 120091) . In the GBM 
bandpass, the light curve is multipeaked with a duration 
of tgo « 70s. GRB 090328 was also detec ted at MeV en - 
ergies by several IPN sat ellites ([G olencts kii et al.ll2009al ). 
including Konus- W^md (Golenetskii et a l.l l2009 d'). 

A spectrum consisting of the first ~ 30 s of GBM data 
(containing the brightest part of the high-energy light 
curve) is well fit by a Band function with a = —0.93 i fc 
0.02, ^ = -2.2±0.1, Ep = 653±45keV (IRau et al.ll2009l) . 
The resulting fluence in the 8keV to 40 MeV (observer 
frame) bandpass is = (9.5 ± 1.0) x lO^^ergcm"^. 
These results are all consistent with the analogous val- 
ues derived by the Koims- Wind instrument, which in- 
clude data out to to + 73 s. Given the wider bandpass 
of the GBM and the higher precision of its fluence mea- 
surement, we adopt this value for remainder of this work. 
Using the observed redshift of 0.736 (§ 12.1.31) . we find a 
prompt fluence of fj — (7.3 ±0.8) x lO^^ergcm^^ in the 
1-10^ keV rest-frame bandpass. 

GRB 090328 was also detected at GeV energies by the 
Fermi-LAT. Much like GRB 090323, flux in the LAT 
band from GRB 09032 8 is detected out to - 900 s af- 
ter the GBM trigger (.Cutini et all l2009f ). Many pho- 
tons with energies above 1 GeV are detected from the 
direction of GRB 090328; however, they ar e all detected 
well after the end of the MeV emission (jCutini et al.l 



120091: IPiron et al.ll2009f ). The highest energy photon de- 
tected from the direction of GRB 090328 was measured 
at ~ 5 GeV {to + 798 s) , while the most energetic photon 
detected du ring the prompt e mission had E « 700 keV 
at to + 60s HPiron et al.ll2009( ). 

2.2.2. Afterglow Observations 

The Swift XRT began target-of-opportunity observa- 
tions of GRB 090328 at 01:26 on 29 March 2009 (- 
15.9 hours after the GBM trigger). A fading X-ray source 
at a = 6h02'"39.''58, S = -41°52'57'.'5 (J2000.0; 2'.'0 con- 
tainment radius) was promptly identified as the X-ray 
afterglow (lKenneall20()9l: iKennea et al.ll2009bD . The re- 
sulting X-ray light curve, showing the afterglow evolution 
over the following 10 days, is plotted in Figure |4l 

The Ultr a violet Optical Telescope (UVOT; 
IRoming et al.l I2005D onboard Swift began settled 
observations of the field of GRB 090328 approximately 
16 hours after the Fermi trigger. The optical afterglow, 
detected in both the [/-band and white filters, was 
identified shortly thereafter (.Oates 2009) . UVOT obser- 
vations of the afterglow continued for ~ 2 weeks, almost 
exclusively in these two filters. We have downloaded the 
UVOT ;7-band data from the HEASARC archive^ and 
conducted photometry with th e images following the 
technique described by iLi et al.l (2006) . Given the large 
uncertainty associated with flux calibration for the white 
filter, we have not included data taken in this band. The 
results of this analysis, not including a correction for 
foreground Galactic extinction {E[B — V] — 0.057 mag; 
iSchiegel et al. 1998), are shown in Table [TUl 

We obtained a single i' image of the afterglow of 
GRB 090328 with GMOS-S on 29 April 2009. This im- 
age was reduced in an identical manner to our observa- 

See http://heasarc.gsfc.nasa.gov. 
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VLA Radio Observations of 
GRB 090328 



Fig. 5. — GMOS-S optical spectrum of the afterglow of GRB 090328. We identify a series of forbidden O emission lines and Fe, Mg, and 
Ca absorption features at a common redshift of 0.7357 it 0.0004. We have only plotted the wavelength range 4500-9000 A, due to the lower 
S/N in the blue (spectrograph throughput) and the red (fringing). 

4000-8000 A), and the second with a central wavelength 
of 8000 A (providing coverage of - 6000-10000 A) The 
spectra were reduced in the manner described in § 12.1.31 
Flux calibration was performed relative to the standard 
star LTT 7379 taken with identical instrumental setups 
on 4 August 2009. 

The resulting spectrum of GRB 090328 is shown in Fig- 
ure [5l Superposed on a relatively flat continuum, we 
identify strong emission lines of oxygen ([O II] A3727, 
[O III] AA4959, 5007) as well as a series of narrow absorp- 
tion features (Fe II A2586, Fe II A2600, Mg II AA2796, 
2803, Mg I A2852, and Ca II AA 3934, 3969), aU at 
z = 0.7357 ±0.0004. 

At this redshift, the prompt isotropic gamma-ray en- 
ergy release from GRB 090328 in the 1-10 keV bandpass 
is E~f,iso = (1.34ib0 . 14) X 10^^ erg. Using the formulation 



Date 


At 


V 




(UT) 


(days) 


(GHz) 


(MJy) 


Mar 30.02 


1.62 


8.46 


172 ± 57 


Mar 30.99 


2.59 


8.46 


337 ± 60 


Apr 1.02 


3.62 


8.46 


783 ± 57 


Apr 3.03 


5.63 


8.46 


195 ± 34 


Apr 4.01 


6.61 


8.46 


674 ± 155 


Apr 5.09 


7.69 


8.46 


214 ± 41 


Apr 6.02 


8.62 


8.46 


523 ± 57 


Apr 6.97 


9.57 


8.46 


809 ± 39 


Apr 11.14 


13.74 


8.46 


603 ± 67 


Apr 14.00 


16.60 


8.46 


643 ± 59 


Apr 16.99 


19.59 


8.46 


886 ± 58 


Apr 24.97 


27.57 


8.46 


410 ± 59 


May 9.96 


42.56 


8.46 


78 ± 43 


May 11.86 


4AA6 


8.46 


21 ± 64 


May 17.92 


50.52 


8.46 


-20 ± 40 



of iLithwick fc S an (2001) , the lower limit on the outflow 
Lorentz factor (assuming a nonthermal spectrum up to 
£;obs « 700 keV) is Fo > 200. 



tions of GRB 090323, and the zero point was calculated 
using predetermined values from the Gemini websiteF^. 
We have also included the simult aneous GROND opti- 
cal and NIR measurements from lUpdike et al.l (|2009al) 
in our modeling, and therefore display them in Table [TU] 
as well. 

Finally, we began observing the aftergl ow of 
GRB 090328 with the VLA on 30 March 2009 (Frai l et al.l 
l2009f l and continued for nearly two months. The data 
were reduced in a manner identical to that described in 
§ 12.1.21 The results of our radio campaign are displayed 
in Table [2] (individual epochs) and plotted in Figure |4] 
(combined epochs). 

2.2.3. Optical Spectroscopy 

We began spectroscopic observations of GRB 090328 
with GMOS-S at 00:0 5 on 30 April 2009 (~ 38.5 hours 
after the GBM trigger: [Ce^ko et al.l [20091) . We obtained 
2 X 1500 s spectra with the R400 grating, the first with 
a central wavelength of 6000 A (providing coverage of ~ 

http : //www . gemini . edu/sciops/iiistruments/ginos/?q=sciops/ 
instrument s/gmos. 



2.3. GRB090902B 
2.3.1. High-Energy Properties 

At 11:05:08.31 on 2 September 2009, the 
Fermi-GBM triggered and located GRB090902B 
(jBissaldi fc Connaughtonll2009D . In the GBM bandpass, 
the light curve consisted of a bright, multipeaked 
pulse with a duration igg « 21s. GRB090902B was 
also detected in a similar bandpass by Suzaku-WAM 
(|Teradaet al.|[2009l ). 

Furthermore, GRB 090902B was bright enough to 
be detected at Ge V energies by the i^ermi-LAT 
(|de Palma et al.ir2009f ) . Like the other events in our sam- 
ple, LAT emission was seen out to -\- 1000 s, at late 
times decayi ng like a power law with temporal index 
aLAT « 1.5 (|Abdo et al.l l2009a). GRB090902B is re- 
sponsible for the highest energy photon detected from a 
GRB to date, with E = 33.4l|^ GeV at io + 82s. Within 
the MeV prompt emission phase, the highest energy pho- 
ton w as measured ai E = 11.16tQ;5g GeV (|Abdo et all 
[2009al) . 

Unlike the other events considered here, the prompt 
high-energy spectrum of GRB 090902B is not adequately 
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Fig. 6.— The broadband radio (blue), optical (red), UV (cyan), and X-ray (black) light curves of GRB090902B. The best-fit model is 
plotted in solid lines (see Table |6] for parameters). The identical model parameters for an isotropic explosion are plotted as the dashed 
lines. The strength of the possible modulation of the radio afterglow caused by interstellar scintillation (e.g.. Frail ct al. 2000a) is indicated 
by the light-blue shaded region. As suggested bv iPandev c t al. (2010), the early {t < O.Sday) optical (and, to a lesser extent, radio) data 
are significantly larger than predicted by our forward-shock afterglow models. The most likely explanation is the presence of reverse-shock 
emission. The lack of a bright radio afterglow at late times suggests the presence of a jet break at t fa 6 days. 



described by a single Band function. While the 100- 
1000 keV bandpass is reasonably well fit by such a model, 
the spectrum exhibits excess emission at b oth low (< 
lOOkeV) and high (> lOMeV) energies. lAbdo et all 
((2009a') have suggested that the high-energy spectrum 
of GRB 090902B can be reproduced as the sum of two 
components: a Band function peaking at ^ 700 keV, 
and a single power-law (photon index T = /3 + 1 — 
1.93) extending over the entire GBM-I-LAT bandpass. 
In this model, the power-law component accounts for 
~ 24% of the total lOkeV to lOGeV fluence. The phys- 
ical mechanism responsible for this complex spectrum is 
still not entirely understood; possible explanations in- 
clude a hadronic origin [eithe r proton synchrotron ra- 
diation ( Razzaoue et alJ 120091) or photohadronic inter- 



actions ( Asanoet al.J 2009|)] or thermal emission from 



the jet phot osphere (jMeszaros fc Ree£ll2000l : lRvgil2004l : 
IRvde et al.l[2010h . 

Using their two-component spectrum, lAbdo et al.l 
(|2009al ) report a total fluence of = (4.59 ± 0.05) x 
10~*ergcm^'^ in the 8keV to 30GeV (observer frame) 
bandpass. At a redshift of 1.8229 (§ [M3D, this 
corresponds to a prompt fluence of (3.83 ± 0.05) x 
ergcm^^ 



10 * erg cm ^ in the 1-10"* keV rest-frame bandpass 



2.3.2. Afterglow Observations 

The Swift XRT began target-of-opportunity observa- 
tions of the field of GRB 090902B at 23:36 on 2 Septem- 
ber 2009 (- 12.5 hr after the GBM trigger). A fading 
X-ray source at a = 17h39™45.''26, 6 = +27°19'28'.'l 
(J2000.0; 2'.'1 containment rad ius) was promptly iden- 
tified as the X-ray afterg l ow (iKennea fc Strattal 120091 : 
lEvansI l2009t iStratta et &L '2009). The resulting X-ray 
light curve is shown in Figure [6l 



The Swift UVOT began concurrently observing the 
field of GRB 090902B and first reported the detection of a 
candid ate optical afterglow con sistent with the X-ray po- 
sition ((S^nso^^^iSiegelHOO^. Subsequent observations 
revealed that the candidate had faded, confirming its as- 
sociation with GRB090902B (Swenson & Stratta 2009). 
In Tabic [Tl] we present UVOT ?7-band observations of 
GRB090902B, reduced in an identical manner to those 
described in § [2X2l 

We obtained a single r' image of the afterglow of 
GRB 090902B with GMOS-N on 3 September 2009. This 
image was reduced in the same manner as those of the 
other events, and the resulting photometry is presented 
in Table [TlJ We have include d optical and NIR p ho- 
tometry of GRB 090902B from iPandev et all ()2010D in 
our modeling, and therefore display them in Table [Tl] as 
well. 

Finally, we began observing the afterglow of 
GRB090902B with th e VLA on 3 September 2009 
(jChandra fc Fraiill2009( ) and continued for 5 months. 
The data were reduced as described in § 12.2.21 The re- 
sults of our radio campaign are displayed in Tables [3] 
(individual epochs) and plotted in Figure [6] (combined 
epochs). We have also included the early radio detection 
at 4.8 GHz reported by the WSRT (jvan der Horst et alJ 
[2009h in Table El 

2.3.3. Optical Spectroscopy 

We began spectroscopic observations of GRB 090902B 
with GMOS-N at 06:29 on 3 September 2009 (^ 
19.4 hours after the GBM trigger: ICucchiara et aO 2009). 
We obtained 2 x 900 s exposures, both with the R400 
grating and a central wavelength of 6000 A, providing 
coverage of ~ 4000- 8000 A. Th e spectra were reduced as 
described in §[2TL3]and ii l2.2.3l Flux calibration was per- 
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Fig. 7. — GMOS-N optical spectrum of the afterglow of GRB 090902B. We identify a series of absorption features from neutral and singly 
ionized Mg, Mn, Fc, and Cr at a common rcdshift of 1.8229 ± 0.0004. We have not plotted the shortest wavelengths due to the decreased 
spectrograph throughput. 



TABLE 3 

Radio Observations of GRB 090902B 



Date 


At 


V 




Facility 


(UT) 


(days) 


(GHz) 


(MJy) 




2009 Sep 3.77 


1.31 


4.8 


111 ± 28 


WSRT 


2009 Sep 3.94 


1.48 


8.46 


141 ± 39 


VLA 


2009 Sep 7.05 


4.59 


8.46 


13 ±31 


VLA 


2009 Sep 8.05 


5.59 


8.46 


130 ± 34 


VLA 


2009 Sep 10.15 


7.69 


8.46 


10 ±32 


VLA 


2009 Sep 11.05 


8.59 


8.46 


80 ±32 


VLA 


2009 Sep 13.10 


10.64 


8.46 


99 ±31 


VLA 


2009 Sep 14.14 


11.60 


8.46 


71 ± 33 


VLA 


2009 Sep 18.04 


15.50 


8.46 


52 ±32 


VLA 


2009 Sep 19.00 


16.46 


8.46 


89 ±36 


VLA 


2009 Sep 25.09 


22.51 


8.46 


26 ±29 


VLA 


2009 Sep 27.08 


24.50 


8.46 


67 ±29 


VLA 


2009 Oct 7.01 


34.43 


8.46 


38 ±28 


VLA 


2009 Oct 9.01 


36.43 


8.46 


66 ±27 


VLA 


2009 Oct 11.81 


39.23 


8.46 


21 ± 31 


VLA 


2009 Nov 6.90 


65.44 


8.46 


9 ±20 


VLA 


2009 Nov 7.77 


66.31 


8.46 


22 ± 19 


VLA 


2009 Nov 9.00 


67.54 


8.46 


48± 19 


VLA 


2009 Nov 14.88 


73.42 


8.46 


31 ± 21 


VLA 


2010 Mar 20.62 


199.16 


8.46 


18± 16 


VLA 



^ Reference: Ivan der Horst eTaLl l(2009l ). 

formed rela tive to spectra of the standard star Feige 34 
(|Okelll990l ) taken with the same instrumental setup as 
on 1 May 2008. 

The resulting spectrum of GRB 090902B is shown in 
Figure [71 Super-imposed on a smooth power-law contin- 
uum, we identify strong absorption features correspond- 
ing to Mg I A2852, Mg II AA2796, 2803, Mn II A2606, 
Fe II A2600, Mn II A2594, Fe II A2586, Mn II A2576, Fe I 
A2484, Fe II A2382, Fe II A2344, Fe II A2260, Fe I A2167, 
Cr II AA2056, 2062, and Mg I A2026, ah at a common 
redshift of 1.8229 ± 0.0004. 

The lack of evidence for Lya absorption down to 
A < 4000 A implies an upper bound on the redshift of 
the GRB host galaxy of z < 2.3. Together with the 
strength of the above features (in particular Fe I A2167), 
we consider it quite likely that the observed system at 



z = 1.8229 derives from the GRB host galaxy. 

At this redshift, the prompt isotropic gamma-ray en- 
ergy release in the 1-10^ keV rest-frame bandpass is 
Ej.iso = (3.20 ±0.04) X 10^"^ erg . Assuming a non thermal 
spectrum up to Eohs ==11 GeV. IAbdo et al.l (p009ail have 
determined a lower limit to the initial outflow Lorentz 
factor of Fo > 1000. 

2.4. GRB 090926 A 

2.4.1. High-Energy Properties 

GRB090926A triggered the GBM on Fer mi at 
04:20:26.99 on 26 September 2009 (|Bissaldil 120091 ) . The 
light curve consists of a single pulse with duration tgo « 
20 s in the GBM bandpass. The prompt emission was 
sufficiently bright to trigger several additional sa tellite 
instruments, ii icluding Suzaku-WKM (iNoda et al"1 [2009). 
Konus- Wind ( iGolenetsk ii et al.' '2009b'), and RT-2 on 
CORONAS-PHOTON (Chakrabarti et al. 2009j. 

Simultaneously fitting both the GBM and LA T (see 
below) data from to to + 21 s, .Bissaldi et ahl (|2OO90 
report that the spectrum is reasonably well modeled by 
a Band function with a = -0.693 ± 0.009, /3 = -2.342 ± 
0.011, and = 268±4keV. The corresponding lOkeV 
to lOGeV fluence (observer frame) is = (2.47±0.03) x 
10"** ergcm"^. These results differ only slightly from 
the values reported fr om other satellites. Kt z — 2.1062 
(jMalesani et al.l[2009l ). the corresponding prompt ffuence 
in the rest-frame 1-10''' keV bandpass is (1.73 ± 0.03) x 
10"'' ergcm"^. 

GRB090926A was also detected by the i^ermz-LAT. 
The emission in the LAT bandpass lasted at least 400 s, 
with possible indications of significant fiux out to ^ 
1000 s after the GBM trigger. Many photons with E > 
1 GeV were detected from the position of GRB 090926A, 
with the highest ene rgy measured of 20 GeV at to + 26 s 
(jUehara et all [20091 ) . 

2.4.2. Afterglow Observations 

The Swift XRT began observing the field of 
GRB090926A at 17:17 on 26 September 2009 (~ 
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Fig. 8.— Early optical and X-ray afterglow of GRB 090926A. 
After an initial decline in the U and V filters, the afterglow experi- 
ences a prominent rebrightening in all filters, peaking at t Ri 1 day. 
Such behavior is difficult to reconcile with standard afterglow the- 
ory, and may suggest a l ate-time (t 3> Atp R^B) injection of energy 
from the central engine l|Rees MeszaroJ[r998 l. 

13 hours after the GBM trigger). A fading X-ray coun- 
terpart at a = 23*^33™36.n8, 6 = -66°19'25'.'9 (J2000.0, 
I'.'S contain ment ra dius) was promptly identified in the 
XRT data (jVetere e t al. 2009; Vcterc' 2009). The XRT 
observed GRB 090926 A for the next 3 weeks, and we plot 
the X-ray light curve in Figures |8] and HI 

Under the control of Sk ynet, four of the 16 in ch diame- 
ter PROMPT telescopes (|Reichart et all 120051 ) at Cerro 
Tclolo Inter- American Observatory (CTIO) observed the 
Fermi-LAT locahzation of GRB090926A beginning 19.0 
hours after the GBM trigger in the V-, R-, and I- 
band filters. Within the Swift-XBT localization, we iden- 
tified an uncatalogued and fadi ng source as the optical 
afterglow of GRB090926A (Ha islip et all l2009bl ial) . In- 
dividual images were automatically reduced using our 
custom, IRAF-based reduction pipeline and then astro- 
metrically aligned and stacked. We subsequently mea- 
sured the afterglow flux with aperture photometry, where 
the inclusion radius was approximately matched to the 
FWHM of the PSF. PRO MPT contin ued to observe the 
field for ten more nights (jHaislip et a l. 2009e,c,d). 

We also obtained 3 epochs of /- and J-band imaging 
of the afterglow of GRB 090926A using the ANDICAM 
(A Novel Dual Imaging CAMera) instrument mounted 
on the 1.3 m telescope at CTICP^. This telescope is op- 
erated as part of the Small and Moderate Aperture Re- 
search Telescope System (SMARTS) consortiunF^. Each 
epoch consisted of 6 individual 360 s /-band observations 
and 30 individual 60 s J-band observations. Between op- 
tical exposures, the telescope was slightly offset and the 
individual J-band exposures were additionally dithered 
via an internal tilting mirror system. Standard data re- 
duction was performed on these images, including cosmic 
ray rejection, overscan bias subtraction, zero subtraction, 
flat fielding and sky subtraction to correct for the NIR 

http : / /www . astronomy . ohio-state . edu/ ANDICAM. 
^® http : //www . astro . yale . edu/ smarts. 
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Fig. 9.— Late-time afterglow and model of GRB090926A. Both 
the X-ray and optical {g'r'i') bandpasses exhibit a steepening de- 
cline at t ^ todays, strongly indicative of a jet break. 



background and the /-band fringing. For each epoch, the 
individual images were then aligned and averaged to pro- 
duce a single frame in each band with summed exposure 
times of 36 minutes in / and 30 minutes in J. 

Relative aperture photometry was performed on the 
SMARTS data, in comparison with a number of nonva- 
riable sources in the field of GRB 090926A. The /-band 
field was photometrically calibrated by comparison (on 
photometric nights ) with Lan dolt standard stars in the 
field of T Phe ( La ndoltl[T99^ . J-band photometric cal- 
ibration was performed using 2MASS (jSkrutskie et al.l 
1006) field stars. 

SMARTS BVRI observations of the field of 
GRB 090926A were also obtained on two photo- 
metric nights a few months after the GRB occurred (16 
and 18 December 2009). For these observations, total 
summed exposure times amounted to 180 s in BRI and 
120 s in V. The absolute photometry of the field was 
again established based on same-night observations of 
the T Phe Landolt standard stars. These observations 
were then used to provide absolute calibration for the 
PROMPT observations of GRB 090926A. 

We obtained additional late-time imaging of the field 
of GRB 090926A on 19 October 2009 with GMOS-S on 
Gemini South. A total of 600 s of exposure time was 
obtained in the Sloan g'-, r'-, and i'-band filters. The 
data were reduced in the manner described in § 12.1.21 
and calibrated in the same way as the PROMPT and 
SMARTS observations. 

The results of our optical and NIR monitoring cam- 
paign of the afterglow of GRB 090926A, uncorrected for 
the modest amount of Galactic extinction [E{B — V) = 
0.024 mag; ISchleeel et aD[T998j . are shown in Table [T^ 
and Figures [5] and [S] 

Finally, the field of GRB 090926A was observed in the 
radio (5.5 GHz) on 1 October 2009 with the Australia 
Telescope Compact Array (ATCA). No source was de- 
tected a t the afterglow lo cation to a 2a limit of < 
1.5mJy (jMoin et al.ll2009l) . 
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2.4.3. Spectroscopy 

iMalesani et al.l (|2009[ ) obtained a spectrum of the af- 
terglow of GRB 090926A with the X-Shooter instrument 
mounted on the 8 m Very Large Telescope UT2. Based on 
the detection of a damped-Lya system and many strong, 
narrow absorption features, these authors derive a red- 
shift of 2.1062 for the host galaxy of GRB090926A. 

At this redshift, the prompt isotropic gamma-ray en- 
ergy release in the 1-10^ keV bandpass is E'^^iso — (1.89± 
0.03) X 10^'* erg. Assuming a nonthermal spectrum up 
to i^obs = 20 GeV, we infer a lower limit for the initial 
Lorentz factor of Tq > 700 (jLithwick fc Sarill200l . 

3. AFTERGLOW MODELING AND RESULTS 

In the standard "fireball" model of GRBs fe.g.. iPiranI 
|2005[) , a compact central engine (likely either a black hole 
or a rapidly spinning proto-neutron star) drives a coUi- 
mated (6* < 10°), uftrarelativistic (Fq > 100) outflow of 
matter and/or radiation. Some dissipative process within 
the outflow (possibly collisionless shocks) gives rise to the 
prompt gamma-ray emission, converting some fraction 
(r]j = Ej^iso/{Ej_iso + -Eke.Iso) of the total relativistic 
energy to high-energy radiation. 

The outflow is ultimately slowed as it sweeps up and 
shock heats the circumburst medium, and synchrotron 
radiation from electrons accelerated at the shock front 
results in the broadband afterglow. The resulting spec- 
trum is well described as a series of broken power laws 
with three characteristic frequencies: Va, the frequency 
below which the radiation is self-absorbed; Vrm the char- 
acteristic frequency of the emitting electrons; and Vc, 
the frequency above which elect rons are able to cool ef- 
flciently through radiation fe.g.. iGranot &: Sarill2002f) . 

The temporal evolution of the afterglow depends on the 
density proflle of the circumburst medium. We consider 
here two possibilities: a constant-density circumburst 
medium [p{r) oc r°], as would be expected in a n environ- 
ment similar to the interstellar medium (ISM: iSari et aD 
I1998D . and a wind-like environment [p{r) oc r~^], as 
would be the case for a massive-star progenitor that 
shed its outer envelope at a constant rate ([Chevalier &: Lil 
I2OOOI) . 

At early times, the afterglow emission appears isotropic 
to distant observers due to the effects of relativistic 
beaming. However, the outflow slows as it sweeps up 
more and more circumburst material. When T(t) w 1/9, 
lateral spreading of the jet becomes important and ob- 
servers will no t ice "missing" em ission from wide angles 
(jRhoadsl [I999I: iSari et al.l I1999D . This hydrodynamic 
transition manifests itself as an achromatic steepening 
in the afterglow light curve. Measuring the time of this 
jet break (tj) allows us to infer the opening angle of the 
outflow. 

In order to ascertain the total relativistic energy out- 
put from a GRB, we require three measurements: (1) 
Ey^iso, the isotropic energy release in the prompt gamma- 
ray emission, which is inferred from the high-energy flu- 
ence and the associated afterglow or host redshift; (2) 
9, the half-opening angle of the beamed emission, in- 
ferred from the detection of a jet break; and (3) Eke, 
the kinetic energy of the blast wave that is powering the 
broadband afterglow, which can be inferred either via af- 
terglow modeling or, more accurately, from late-time ra- 



dio calorimetry in the nonrelativistic phase (iBerger et al.l 
[2OOI iFrafl et al.l [20051: Ivan der Horst et al.l 120081) . We 

stress here that we are neglecting contributions from non- 
electromagnetic phenomena (neutrinos, gravity waves, 
etc.) and slower-moving material (i.e., supernova emis- 
sion) , and so are providing only lower limits on the total 
energy budget. 

Unfortunately, none of the radio afterglows in our 
sample of LAT events are sufficiently bright to perform 
calorimetry in the nonrelativistic phase. Instead, we con- 
struct afterglow models (including both the standard for- 
mulation and corre ctions for radiativ e losses and inverse- 
Compton emission; ISari fc Esinl[200l and compare these 
with our o bservations using the multi-parameter fitting 
program of lYost et al.l ([20031). Our objective is to trans- 
late the observed three critical frequencies, together with 
the peak fiux density, -Fl/^max, and the jet break time, 
tj, into a physical description of the outfiow. In par- 
ticular, we shall attempt to estimate seven parameters: 
Eke, the kinetic energy of the blast wave; n, the par- 
ticle density of the circumburst medium (or, alterna- 
tively for a wind- like circumburst medium, A^, where 
p = 5 X lO^^A+r"^ gcm^^jHH; fe, the fraction of the to- 
tal energy apportioned to electrons; e^, the fraction of 
the total energy apportioned to the magnetic field; p, 
the electron power-law index; Ay, the host-galaxy ex- 
tinctiorF^: and 9, the jet half-opening angle. 

Optical magnitudes have been converted to flux den- 
sities after correcting for Galacti c extinction using zero 
points from lFukugita et al.l ([19951 ). To account for differ- 
ences in instrumental configurations, we have applied a 
7% cross-calibration uncertainty to all data points before 
calculating the models. All reported uncertainties have 
been determined using a Monte-Carlo bootstrap analysis 
with 1000 trials and represent only statistical errors as- 
sociated with the fit. Systematic errors associated with 
model uncertainties are potentially larger and difficult to 
estimate. 

3.1. GRB 090323 
3.1.1. Preliminary Considerations 

Before proceeding to a detailed model, we derive some 
initial constraints by looking at the spectral and tem- 
poral behavior of the afterglow. Considering first the 
X-ray afterglow, we fit the flux density to a power- 
law decay and find a best-fit index of ax = 1-5 ± 0.1 
(x^ = 28.2 for 22 degrees of freedom, d.o.f.). Perform- 
ing an analogous power-law fit to the X-ray spectrum, 
we find (3x = IMtait ix^ = 8.0 for 13 d.o.f.; see the 
on-line compilation of N.R.B. for details). 

We perform a similar fit in the optical (r', R, and i' 
filters), forcing the decay index to be identical in all fil- 
ters and ignoring the points at i > 10 days due to host 
galaxy contamination (see below) . Though the fit quality 

For a wind-like medium, the circumburst density is normalized 
for a progenitor mass-loss rate of M = 10~^ Mq yr~^ and a wind 
speed of = 1000 kms^^ l lChevalier fc LillTgggh . The units of 
At are thus g cm~^, and the conversion to a particle density (as a 
function of radius) is given by n = 30A»r]~-, cm""^. 

We have assumed an SMC-like extinction curve for all events 
here {Pel 1992; Kann et al. 2006). Given the relatively modest 
amounts of dust inferred for these host galaxies, this choice does 
not affect any of our primary conclusions. 
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TABLE 4 

GRB 090323 Afterglow Best-Fit Parameters 



^'KE.iso 
(10^2 erg) 


A. 
(g cm-l) 


(%) 


{%) 


61 

{°) 


p 


AvChost) 
(mag) 


Xi (d.o.f.) 


2001^° 




6.8 ±0.5 


0.31 ±0.12 


^■"-0.1 


2.81 ±0.06 


0.12 ±0.03 


1.39 (70) 



is much worse (x^ = 178.4 for 34 d.o.f.), we find a sim- 
ilar decay index as seen in the X-rays: ao — 1-8 ± 0.1. 
Likewise, a spectral fit of the multi-color GROND data 
obtained at t = 1.12 days results in a spectral index of 
/3o = 0.95 ±0.01. 

The comparable spectral and temporal indices in the 
optical and X-rays suggests that both bandpasses fall in 
the same synchrotron spectral regime. In other words, 
the cooling frequency should fall either below the optical 
or above the X-ray bandpass for the duration of these 
observations {t w 1-10 days). Examining the "closure" 
relations between a and l3 in different circumburst me- 
dia and spectral regimes (e.g.. [Price et al. 2002) . we find 
that we can rule out a low cooling frequency at large sig- 
nificance, as this would require a = (3/3 — l)/2 1.0 for 
both a wind-like and ISM-like circumburst medium. 

The relatively flat radio light curve (modulo effects 
from interstellar scintillation) until t « 20 days has two 
important implications. First, it suggests a wind-like cir- 
cumburst medium, as the flux would be expected to rise 
in proportion to t^/^ for i^a < v < Vm in a constant- 
density environment. From this and our previous X- 
ray and optical spectral and temporal decay indices, we 
can infer that the electron spectral index is relatively 
steep, V fti 2.7 (e.g iShen et aljr2006l: iStarling et alll2008l: 
iCurran et al.ll2009( ) . Second, together with the lack of ev- 
idence for late-time steepening in the X-ray and optical 
light curves, this suggests that if a jet break is required, it 
should occur at t > 20 days. (Another possibility for the 
decay could be the peak frequency v„i passing through 
the radio bands.) 

Finally, the relatively bright r'-band detection at i = 
130 days is almost certainly dominated by flux from the 
host galaxy and not the afterglow of GRB 090323. If we 
assume this flux is due entirely to host light, the host 
will contribute a significant fraction of the flux in some 
of the late-time points of our optical light curve (~ 40% 
at t = 14 days in r'). We assume a flat spectrum and 
subtract a host contribution of = 0.40 /iJy from all 
optical data points for our afterglow modeling. 

3.1.2. Modeling Results 

With the above constraints in hand, we have modeled 
the afterglow of GRB 090323 with the software described 
above. The resulting best-fit wind model is plotted in 
Figure and the derived parameters are provided in Ta- 
ble 131 The overall fit quality is reasonable (x^ = 97.4 
for 70 d.o.f.), with comparable residuals in the radio, op- 
tical, and X-ray bandpasses. We were unable to find 
any physically plausible models for a constant-density 
medium with reasonable fit quality. 

As expected, the cooling frequency lies at or above 
the X-ray bandpass for the duration of our observa- 
tions {t > Iday). As a result, the ratio of i?KE,iso 
to Aif is required to be larger than inferred for previ- 



ousGRBs (e.g.. iPanaitescu fc Kumaijl2001bl : IVost et afl 
|2003( ). While the isotropic blast -wave kinetic energy 
is relatively large compared to that of previously mod- 
eled GRBs (-EKE.iso = 2.0 X 10^'' erg), it is in fact com- 
parable to the prompt gamma-ray energy release, as 
would be expected for reasonable values of the gamma- 
ray efficiency rj^. The inferred density is slightly lower 
than usual {A^ = 0.12 g cm~^), though smaller values 
have been reported i n the literature (e.g., GRB 020405; 
IChevalier et"aI1[200l . 

The models result in a jet break time of tj = 
17.8l|^^7^ days. This occurs after the X-ray and optical 
observations have stopped (except for the host detection 
in the optical), and is therefore most directly constrained 
by the radio light curve. In Figure[2]we also plot the same 
model parameters but for an isotropic explosion, where 
the turnover in the radio at late times is due to the peak 
frequency Um passing through the radio bandpass. The 
current radio data provide only weak constraints on the 
jet break time upper bound, although the models sug- 
gest a much smaller uncertainty in the opening angle. 
We return to the robustness of our determination of 9 in 

§1331 

While the jet break occur s rather late in the observer 
frame (e.g.. lZeh et al.|[2006[) . the inferred opening angle 
{9 = 2.6°) is still relatively small compared with pre- 
vious samples. Though the dependence is not strong 

{9 oc _(£^KE,isoM*)"^^^), the large ratio of SKE.iso to 
effectively lowers the opening angle for a given jet break 
time. 

After applying the coUimation correction, we find that 
the true energy release of GRB 090323 is = 4.ll^-9 x 

10^^ erg, Eke = '^■^-o.t ^ 10^^ erg. We com pare these 
results with a larger sample of events in § 14.11 

3.2. GRB 090328 
3.2.1. Preliminary Considerations 

Following the results of § 13.1.11 we first consider the 
spectral and temporal indices of GRB 090328 in the ob- 
served bandpasses independently. The X-ray afterglow 
light curve is well fit by a single power-law decay with 
index ax = 1.65 ± 0.09 (x^ = 24.1 for 26 d.o.f.). The 
spectral index, on the other hand, is less well constrained: 
I3x = l.lloi (X^ = 6.7 for 14 d.o.f.). We note that the 
X-ray spectral fit indicates the presence of an additional 
absorption component in the host galaxy of GRB 090328 
{Nh = 4.9l^J X 10^1 cm-2 at z = 0.736). 

The relative paucity of optical data, particularly at late 
times, makes a similar analysis challenging. Fitting the 
UVOT J7-band light curve to a single power-law decay 
results in an index of ao = 1-2 ± 0.1; however, the qual- 
ity of the fit is extremely poor (x'^ = 24.5 for 5 d.o.f.). 
An initially steep decay in the light curve becomes more 
shallow at later times, much like is seen in the UVOT 
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TABLE 5 

GRB 090328 Afterglow Best-Fit Parameters 



-BKE.iso A* ee EB p Ay(host) xl (d.o.f.) 
(1052 erg) (gcm-i) (%) (%) (°) (mag) 

n+l 0.n+l°l 9±2 O.9I02 5. 2^0 7 2.58^0.^7 0.41 ± 0.04 1.63 (43) 



white filter (jMarshall et al.l 120091 ) . This result suggests 
the presence of host-galaxy contamination at late times. 

The best-fit spectral index assuming a power-law 
model for the simultaneous multi-color GROND data is 
(3o = 1.24 ± 0.04, although the fit quality is extremely 
poor (x^ = 85.5 for 5 d.o.f.). The poor fit statistic is 
dominated almost entirely by a single filter (J band), 
however, and the relatively steep spectral index (com- 
pared to Pox ~ 0.8) appears to be robust. This is con- 
sistent with the presence of dust indicated by the X-ray 
spectral fits. 

In the case of GRB 090328, examining the a-/3 clo- 
sure relations yields little insight. Because of the poor 
fit quality in the optical and the large uncertainty in the 
X-ray spectral index, we are unable to rule out much of 
the available parameter space. We note, however, that 
with the exception of the excess absorption, the X-ray af- 
terglow light curve and spectrum of GRB 090328 closely 
resemble those of GRB 090323. Similarly, if we assume 
that the late-time optical data are dominated by host- 
galaxy light, the optical light curve is also compatible 
with that seen from GRB 090323. 

3.2.2. Modeling Results 

The results of our best-fit model for the afterglow of 
GRB 090328 are shown in Figure 2] and the parameters 
are provided in Table [5l The overall fit quality is rea- 
sonable (x^ — 70.1 for 43 d.o.f.). Only the early radio 
light curve, which suffers from strong scintillation, differs 
systematically from the model predictions. 

Much like in the case of GRB 090323, we find that only 
a wind-like circumburst medium can reproduce the ob- 
served light curves with plausible physical parameters. 
While the slow rise in the early radio afterglow is better 
explained by a constant-density environment, our mod- 
eling indicates ISM-like solutions require eg — 1 and 
es — (i.e., a fully radiative solution). 

We note that our solutions require the presence of a 
relatively bright host galaxy — if the late-time [/-band 
points are not predominantly due to host light, the fit 
quality of our models decreases somewhat {xf, 2.0). 

The steep decline in the radio light curve at late times 
offers relatively tight constraints on the presence of a 
jet break. In fact, the break is consistent with (though 
not formally required by) the final few points in the X- 
ray light curve as well. For our best-fit model, we in- 
fer tj = 6.4tJ2g°days and 6 = (5.2^^;^)°. Primarily be- 
cause of the relatively sparse optical data, we find several 
unique models with similar though slightly worse overall 
fit statistics. We note, however, that all these solutions 
result in comparable opening angles, the most important 
parameter for our energetics calculations. 

Correcting for the effects of beaming, we find an energy 
release of = 5.5tti x 10^° erg and Eke = 'i-^tto x 
10^° erg. The total relativistic energy release, Etot ~ 



]^q51 gj.g^ somewhat less than the value we derive for 
GRB 090323. 

3.3. GRB090902B 
3.3.1. Preliminary Considerations 

A single power-law model with ax — 1.36 ± 0.03 pro- 
vides a good fit to the X-ray light curve of GRB 090902B 
over the entire span of the XRT observations (x^ = 94.8 
for 107 d.o.f.). Likewise, the X-ray spectrum is well fit 
by a power law with (3x = 0.90 ± 0.13 (x^ = 42.8 for 56 
d.o.f.). The X-ray spectrum strongly favors the presence 
of absorption in excess of the Galactic value along the 
GRB line of sight. Assuming the dust arises in the GRB 
host galaxy at z = 1.82, the best-fit host column density 
is Nh = (6.3^1:7) X 10^1 cm-2. 

The decline inferred from the early-time ROTSE-III R- 
band observations is significantly steeper than the late- 
time ( t > 1 day) optical decay. Following iPandev et al.l 
(l2oTnl) . we attribute this emission largely to the presence 
of a reverse shock and do not include these data points 
in any subsequent analysis. Simultaneously fitting the 
remaining optical data to a power law of the form f^, oc 
t-'^v-P, we find ao = 0.89 ± 0.05 and Po ^ 0.76 ± 0.07 
(X^ = 36.4 for 18 d.o.f.). 

The large discrepancy between the X-ray and opti- 
cal temporal decay slopes suggests that the two band- 
passes fall in different synchrotron spectral regimes (i.e., 
1^0 < i^c < t^x)- Furthermore, the relatively shallow 
optical decline strongly rules out a wind-like circum- 
burst medium. We find that both the X-ray and op- 
tical bandpasses are broadly consistent with expansion 
into a constant-density medium and an electron index 
p w 2.2. The observed X-ray temporal decline is some- 
what steeper than what is predicted, but this could be ac- 
counted for by synchrotron radiative losses at later times. 

The radio light curve appears to decay from the earliest 
observations at to + 1.5 days onward. This behavior is not 
expected in the simple forward-shock model where radio 
afterglows are typically observed to rise (for a constant- 
density medium) or remain flat (for a stratified environ- 
ment) until either a jet break occurs or the peak fre- 
quency reaches the radio bandpass (typically weeks after 
the burst). 

However, bright, early-time radio emission is actually 
rather common, having bee n detected in a quarter of 
all radio afterglows to dat e (jSoderberg fc Ramirez-Rui^ 
120031 : [Chandra et al.ll2010l ). This prompt component is 
most commonly attributed to the reverse shock propa- 
gating in the shocked ejecta. The reverse-shock emis- 
sion fades away quite quickly, and at late times {t > 
5 days) the afterglow will be dominated by the forward- 
shock emission. GRB 990123 was the first known ex- 
ample (jKulkarni et al.|[T99 9l. but t here have been many 
more claimed detection s since (e.g., jBerger et aLll2003bl: 
iPanaitescu.2005; .Chandra et al.ll2010[ ). 
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TABLE 6 

GRB 090902B Afterglow Best-Fit Parameters 



-BKE.iso " 
(10^2 gjg) (cm-^) 


(%) 


(%) 


e 

n 


p y4t/(host) 
(mag) 


(d.o.f.) 




15tt 




o -+0.4 
■^■^-0.3 


r, 99 + 0.08 r, 1O+0.08 

^•^^-0.04 "--^"-O.OS 


1.06 (137) 



The radio light curve of GRB 090902B most closely re- 
sembles that of GRB 991216 in which the early power-law 
decline is due to the o verlap of emission from the reverse 
and forward shocks (jFrail et al.l l2000bl) . Since we are 
concerned here with the modeling of the behavior of the 
forward shock, we will exclude the early radio point and 
postpone a full discussion to a subsequent paper dealing 
with a co mpilation all known GRBs with prompt radio 
emission (IChandra et al.ll20Tol) . 

3.3.2. Modeling Results 

Using the constraints derived above, our best-fit model 
for GRB 090902B assuming a constant-density circum- 
burst medium is shown in Figure |6l and the derived 
model parameters are presented in Table [51 The over- 
all fit quality is quite reasonable (x^ = 145.8 for 137 
d.o.f.), although the fit quality is somewhat better in 
the X-ray and radio than in the optical/UV. We stress, 
however, that we have not included the first two i?-band 
data points and the first VLA detection in our fitting, 
as we believe the flux at these points is dominated by 
reverse-shock emission. 

Much like the other events considered here, our inferred 
parameters for GRB 090902B suggest a large afterglow 
kinetic energy (-Eke.iso = G.St.Q Q x 10^^ erg; a factor of 
flve less than the prompt gamma-ray energy release) and 
a low circumburst density (n = 5.8li[ g x 10~'*cm~^). 
In fact, we are unable to find any acceptable solutions 
(xt < 2) with n > 10~^ cm~^. The inferred circumburst 
density for many long-duration GRBs is lower than typ- 
ical values observed in dense molecular clouds (where 
presumably their massive-star progenitors have formed). 
Yet the value we have derived for GRB 090902B is closer 
in fact to what one might find in the ambient ISM or even 
the intergalactic medium (IGM), as has been found for 
the f ew short-hard GRBs with sufficient af terglow data 
(e.g.. iPanaitescul [20061 : IPerlev et al.ll2009bf ). We return 
to this issue in more detail in § 14.31 

Neither the X-ray nor the (less constraining) optical 
bandpass formally require a jet break over the duration 
of our observations. In this case, like GRB 090323, we 
must rely largely on the radio. Our models suggest that 
the peak synchrotron frequency, Vrm will not fall to the 
radio bandpass until t > 40 days. The declining radio 
light curve at late times therefore requires a jet break 
at tj — 6.2lQ g days. Again, because of the large ratio 
of -EKE.iso to n, this corresponds to a relatively narrow 

opening angle: 6 = {3Atgi)°. 

Correcting for coUimation, we find the true energy re- 
lease from GRB 090902B to be = (5.6±1.5) x 10^^ erg, 
£;ke = 1.21;^:Ix 10^1 erg. 

3.4. GRB090926A 
3.4.1. Preliminary Considerations 



Fitting a single power-law decay to the X-ray light 
curve, we find a reasonable quality fit (x^ — 120.8 for 97 
d.o.f.) with ax = 1-43 ± 0.03. The fit quality improves 
somewhat if we allow for a steeper decay at late times, 
although the post-break decay index is not particularly 
well constrained: ax,i = 1-35 ± 0.04, ax,2 = 
tb = 9.1^1-4 days, = 111-9 (95 d.o.f.). The X-ray spec- 
trum is well fit (x^ = 41.7 for 51 d.o.f.) by a single power 
law with index /3x = 1.12 ± 0.13 and does not require 
any absorption in addition to the Galactic component. 

The behavior in the optical bandpass, however, is sig- 
nificantly more complex (Fig. [8]). The flux in the ini- 
tial UVOT U- and X^-band observations declines until 
t « 0.8 day, at which point a rebrightening or plateau 
phase is evident in all five observed filters (UBVRI). 
The optical light curve of GRB 090926A is reminiscent 
of the rebrightening see n at t « 1 day from GRB 970508 
(iDiorgovs ki eTall [l997l ) and GRB 071003 (jPerlev et al.l 
2008 and references therein). Such a phenomenon is dif- 
ficult to reconcile with the standard afterglow paradigm. 
Possible explanations in clude a sharp chang e in the cir- 
cumb urst de nsity (e.g.. iLazzati et ahl l2002t iTam et all 
120051 but, c.f. iNakar fc Granotil2007i ror a smooth injec- 
tion of energ y into the forward sho ck from the central 
engine (e.g., iRees fc MeszarosI 119981 ). A more detailed 
analysis of the early-time behavior of GRB 090926A is 
beyond the scope of this work. 

The optical light curve appears to peak at i ss 1 day 
in aU bands. Fitting the PROMPT BVRI data at 
t > I day to a single power law results in a best-fit index 
of ao = 1.38 ±0.02 (x^ = 181.6 for 119 d.o.f.). However, 
extrapolating these results to the late-time Gemini imag- 
ing at t w 23 days greatly overestimates the observed fiux 
(Fig.|9]). Much like the X-ray data, this strongly suggests 
a steepening of the optical light curve at i > 8 d. 

We have also performed a joint spectral and temporal 
fit to the PROMPT data to estimate the optical spectral 
index, /Sq- We include only the V-, R-, and /-band data, 
as both the U and B bands are likely affected by Lya 
absorption at z = 2.11. We then find a best-fit optical 
spectral index of (3o — 1.03±0.05. The formal fit quality 
is relatively poor, however (x^ = 183.2 for 68 d.o.f.). 
Finally, the optical to X-ray spectral index at t > 1 day 
is l3ox ~ l-l- 

Taken together, we find similar temporal {ax,i = 1-35, 
ao,i = 1.38) and spectral {f3x = 1-12, f3o = 1-03) indices 
in the X-ray and optical for t « 1-9 days. Furthermore, 
the optical to X-ray spectral index, Pox, is comparable 
to both the optical and X-ray spectral indices. These 
facts strongly suggest that both the X-ray and optical 
bandpasses fall in the same synchrotron spectral regime. 

Considering the various synchrotron closure relations, 
the best fit appears to occur when both bandpasses fall 
above the cooling frequency, Vc- The afterglow decay in 
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TABLE 7 

GRB 090926A Afterglow Best-Fit Parameters 



-EkE.Iso 
(10^2 erg) 


A* 
(gcm~i) 


(%) 


(%) 




n 


p Av{host) 
(mag) 


xl (d.o.f.) 


24.2tl'i 


-0.02 


301^ 


n 07+0. 07 


7.1 ±0.2 


2.19 ±0.02 0.06 ±0.01 


1.18 (149) 



this regime is independent of circumburst medium. Such 
a low cooling frequency is somewhat unusual, though not 
unprecedented (e.g., GRB 050904; iFrail et all |2006[) in 
GRB afterglows. The implied electron index in this case 
would be p « 2.3. 

Finally, we note that the steepening in both the X-ray 
and optical bands at i « 9 days is strongly suggestive of 
a jet break. We examine this possibility in greater detail 
in the following section. 

3.4.2. Modeling Results 

Unfortunately, without a radio light curve we cannot 
uniquely solve for the physical parameters of the after- 
glow of GRB 090926A. We can, however, at the very least 
more robustly constrain the jet break time and opening 
angle of this event, as well as identify broad trends in 
the parameters associated with acceptable solutions. As 
discussed previously, we consider only data at t > 1 day 
due to the rebrightening observed in the optical before 
then. 

As expected, we were able to find a number of accept- 
able solutions, both for constant-density and wind-like 
environments. One example for a wind-like circumburst 
environment is shown in Table[7]and FigurelH The model 
provides a reasonable description of the data in all band- 
passes, with = 175.1 (149 d.o.f.). Without radio data, 
however, the solution presented in Table [7] is not unique 
in a global sense, and therefore should only be taken as 
representative of a larger family of solutions. 

Most importantly for our purposes, a jet break at 
t « 9±2 days is required to provide a reasonable fit to the 
late-time X-ray and (particularly) optical data (Fig. |9]). 
The precise translation into an opening angle is depen- 
dent on the circumburst density (both profile and abso- 
lute value) and the afterglow energy. Taking the values 
derived above into account, we find a typical value for 
a wind-like circumburst medium oi 9 w 7°. Based on 
the observed spread of -Eke and A, in our wind mod- 
els, we adopt an approximate uncertainty in this value 
of 6* R:! (7tf)°. 

Correcting for coUimation, we find a prompt energy 
release of = 1.4:~^q'^ x lO'^^erg. Assuming a reason- 
able value for the gamma-ray efficiency (77-, > 10%), we 
adopt a lower limit on the isotropic blast- wave energy of 
-EKE.iso ^ 10^^ erg. We then find a coUimation-corrected 
afterglow energy of Eke 5 x 10^° erg. 

4. DISCUSSION 

4.1. Central- Engine Constraints I: Energetics and 
Remnants 

In Table |51 we summarize the primary results from 
this work, including the redshift, initial Lorentz factor, 
beaming angle, density, and coUimation-corrected energy 
release for each of the four LAT GRBs considered here. 



Before the launch of the Swift satellite in 2004, 
the majority of all well-observed afterglows were 
inferred to be highly co l limate d, with opening an- 
gles e < 10° (jZeh et all |2Q06f). The most notable 
exceptions were the handful of the most nearb; 
events, includ i ng G RB 980425 (Galama et al 199 
Kulk arni et all [l998l) GRB 03 1203 (ISoderberg et al l 
[2OOI 



2004; 



2004 : iSazonov et all 120041: iT homs en et a" 
Cob b et al.l 120041: iMalesani et al .^ 2004: Gal- Yam et al 
[2004), a nd (yost-Swi ft ) GR B 060218 (iMirabal et al 

m 



Ca mpana et al.l 



2006t 



200a 



ISoderberg et aLl 

iModiaz et al.' '2006; 'Sollerman et al.' '2006': ' Plan et al] 
12006: Fcrrcro et al, 2006; Maeda et al.. 200% all of 
which appear to be isotropic explosions that were 
energetically dominated by their nonrelativistic ejecta 
(i.e., their associated supernovae). 

The typical GRB afterglows discovered by Swift, how- 
ever, did not fit neatly into this simple bimodal pic- 
ture. The afterglows of most Swift GRBs, including the 
large fraction at cosmological distan ces, exhibit a much 
broad er range of openin g angles (jKocevski fc Butleil 
l2008t iRacusin et a"Lll2009l) . with some extreme events 
lacking a detectable jet break signature in the X-rays out 
to hundreds of days aft er the burst (e.g., GRB 060729; 
IGrupe et al1l2?M l20Tfll ). 

All four of the LAT-detected events we have studied 
here are consistent with a relatively high degree of coUi- 
mation {6 < 10°). In this respect, then, the afterglows of 
LAT-detected GRBs more closely resemble the pie-Swift 
sample. This in and of itself is not entirely surprising, 
given the tremendous high-i?..y iso bias for events detected 
by the LAT (Figured]). 

In Figure 1101 plot the two-dimensional (prompt 
-|- afterglow) coUimation-corrected relativistic energy re- 
lease for the four LAT events in this work, compared with 
several additional samples. Shown in red are 11 pie- Swift 
GRBs at cosmological distances {z > 0.5) for which suf- 
ficient broadband (X-ray, optical, and radio) data exist 
to constrain both E.y and Eke- We find the logarithmic 
mean for the sum of these two values, i^rei ~ E^-\-EkeF^ 
to be (i?roi) = 2.0 X 10^^ erg, with an error in the mean of 
0.17 dex. The solid black line indicates a constant value 
of Erci corresponding to this mean, while the gray shaded 
regions indicate Icr, 2a, and Scr confidence intervals. 

The gray dotted line in Figure [10] represents a con- 
stant gamma-ray efficiency of 77^ = 50%. Most events 
are roughly consistent (or even exceed) this value. This 
result presents a problem for most internal shock mod- 
els of the prompt emission, which pr edict a maximal 
gamma-ray efficiency of y-y < 1 0% (jKobavashi et al.l 
ll997HDaigne fc Mochkovitchll 19981 ). Alternatively, if the 

23 Assuming an approximately log-normal distribution, we calcu- 
late the weighted mean of (logjg(i?tot )), where a = SEtot / (In(lO) X 
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TABLE 8 

Summary of GRB Parameters 



GRB 


z 


(10 4 erg cm 2) 


fp a 

{10^^ erg) 


To 


e 

n 


(10^1 erg) 


Eke 
(lO^i erg) 


n/A, 
(cm^-^ / gcm-i) 


090323 
090328 
090902B 
090926A 


3.568 
0.7357 
1.8229 
2.1062 


1.50 ±0.20 
0.73 ± 0.08 
3.83 ± 0.05 
1.73 ± 0.03 


3.99 ± 0.53 
0.134 ±0.014 
3.20 ±0.04 
1.89 ± 0.03 


> 600 

> 200 

> 1000 

> 700 


2 6+''''^ 

5 2+l'4 
^•^-0.7 

04+0.4 

7+3 
'-1 


4 1+2.9 
^■^-0.8 
5^+0 '13 
"■°'^-0.18 

5.6 ± 1.5 

14tf 


2 1+2-4 
^■^-0.5 

1-^-0.4 

^ 0.5 


12+°-02 

11+0-03 
^■11-0.05 

(5.8j:l-j) X 10-4 

. . . b 



^ 1—10 keV observer frame bandpass. 

^ We have not included an estimate for the density of GRB 090926A, as this parameter was not well constrained due to the 
lack of radio data. 



early afterglow undergoes a relatively long-lived radia- 
tive phase, we may be si gnificantly underestim ating the 
blastwave kinetic energy (IGhisellini et al.llMol) . 

As discussed previously, the three most nearby GRBs 
are clearly subenergetic outliers. The energy release 
from the supe rnovae associated with these events, fisN ~ 
1051-10^2 gj.g (jiwamoto et al.lUOQSl: [Mazzah et al.|[2006l ). 
is several orders of magnitude larger than the energy ap- 
portioned to the relativistic ejecta. To further distinguish 
these events, the rate of such subluminous GRBs suggests 
they are many times more common (p er unit volume) 
than the typical cosinological GRBs ( Soderberg et alj 
200a ICobb et al.l[2006t iGuetta fc Vallell2007D . although 



Butler et al. I (|201(t r have suggested that both samples 
may be described by a single luminosity function — that 
is, subluminous GRBs may simply extend continuously 
from the higher-energy population. 

Of the LAT events studied here, the total relativistic 
energy output from only a single one (GRB 090328) is 
consistent with the pre- Swift mean at the 3cr level. This 
is not entirely surprising, however, as the E^^i distribu- 
tion of pre- Swift GRBs exhibits a reasonable dispersion 
(^ 0.55 dex, or a factor of 3.5). It is clear, however, that 
GRB 090323, GRB090902B, and GRB090926A aU fall 
at the very high end of the pre-Swift distribution. In 
particular, the only event comparable to GRB 090926A, 
with « 10^^ erg and Eke relatively unconstrained, 
is GRB 970508 (S^i ~ 1.5 x 10^^ er g, dominated by 
the large afterglo w kinetic energy; lYost et all 120031 : 
iBerger et al1[200l . 

GRBs 090323, 090902B, and 090926A appear consis- 
tent instead with a subs a mple of the brightest Swift 
events from iCenko et ahl ()2010l ) (see also Frai l et al.l 
[200aiChandra et al.ll2008l ). This, too, is not unexpected, 
as these Swift events were chosen on the basis of large 
^'7,iso values, and therefore should in large part mimic at 
least this component of the LAT selection effects. 

Much like several events in the bright Swift sample, 
GRB 090926A appears to exceed the canonical GRB rel- 
ativistic energy release of 10^^ erg by roughly an order 
of magnitude. We stress that our methodology provides 
a relatively conservative estimate for the prompt energy 
release of GRB090926A, for a number of reasons. First, 
we consider only the rest-frame 1-10* keV bandpass for 
all events considered in Figure 1101 Extrapolating mea- 
surements from previous instruments, with bandpasses 
typically extending only to ~ 1 MeV (only a few hundred 
keV for Swift) , up to the GeV range would introduce sig- 
nificant uncertainties in Ey^iso- While allowing for a more 
robust burst-to-burst comparison, we have not included 



a significant fraction of the detected high-energy fluence 
in our energy calculations for these LAT events (~ 30% 
for GRB090926A). 

Second, wc have derived the opening angles of all 
our events using parameters inferred from our broad- 
band modeling. Consequently, the i^rei we derive for 
GRB 090902B is a factor of several lower than that re- 
ported by other authors (§ 14.31 and 14. 4p . Most im- 
portantly, for GRB 090926A we have deliberately used 
only wind-like models to derive the beaming angle, even 
though our data are insufficient to rule out a constant- 
density medium. For ISM-like models, the opening an- 
gles we found were on average 75% larger, corre- 
sponding to a beaming-corrected prompt energy release 
of sa 4 X 10^2 gi-g. 

Together with the discovery of several comparable 
events to GRB 090926A in the past few years, includ - 
ing GRB 050904 {E.^i « 2 x 10^^ erg; EiiLeL al.lt2006|), 
GRB 070125 {E.^i « 3 x 10^^ erg: IChandra et al. 2008), 
GRB050820A {E,^i w 4 x 1052gj.g. |Cenk o et al. 2010), 
and GRB 090423 {E.^i > 1 x lO^^erg; Chandra et all 
l2oToh . we now believe there is substantial evidence in 
favor of a subpopulation of GRBs with relativistic en- 
ergy outputs either very near or above 10^^ erg. We refer 
to such events as hyper-energetic GRBs in what follows, 
and outline some of the implications of this particular 
energy threshold. 

Much as was arg ued by iStarhng et aU (|2009D for the 
case of GRB 080721, the total energy budget is an im- 
portant diagnostic for any central-engine model. In par- 
ticular, models for which the outflow is powered by the 
spin-down of a highly magnetized {B > 10^^ G) proto- 
neutron star are subject to strict constraints on the total 
energy budget of Stotai < 3 x 10^^ erg (the rotational en- 
ergy of a maximally spinning 1.4 Mfr, neut ron star; e.g., 
[Thompson et al.l 120041 : iMetzger et al.ll2007| ). 

Only accounting for the relativistic energy output, 
GRB 090926A seems to approach within at least a factor 
of a few of this limit. While we caution that there are still 
significant uncertainties associated with the models used 
to infer the afterglow parameters (and hence i^roi), we 
have not yet accounted for additional sources of energy, 
including (nonrelativistic) SN emission, radiative losses 
at ear l y times due to brig ht X-ray flares (jBurrows et al.l 
l2005at iFalcone et al.ll2007i ). and synchrotron losses dur- 
ing t he later phases of afterglow evolution ()Yost et al.l 
mm- Even if we have overestimated the relativistic en- 
ergy output for these events by a factor of several and 
the magnetar energy limit is not strictly violated, the 
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tremendous efficiency required by this process strains 
credulity. 

4.2. Central- Engine Constraints II: Lorentz Factor, 
Opening Angle, and Acceleration Mechanisms 

Unlike magnetar models that are powered by the spin- 
down of a highly magnetized proto-neutron star, GRB 
models in which the core of an evolved massive star 
collapses promptly to form a black hole and accretion 
disk syste m ("Type I collapsars" following the nomen- 
clature of iMacFadven et al.ll200l1 ) have significantly re- 
laxed constraints on the total energy budget that are 
in at least some cases capable of accommodating hyper- 
energetic events. For instan ce, the initial models of 
IMacFadven fc WooslevI (|1999( ) that begin with a 35 Mq 
He star with a 10 Mq evolved core lead to accretion rates 
of ~ 0.1 M© s~^ onto a 3 M© rotating black hole. Assum- 
ing approximately continuous feeding over the fallback 
time of the stellar envelope (~ 10s), this corresponds to 
a total accreted mass ~ 1 M© . In addition to the accre- 
tion process, an even larger reserve lies in the rotational 
energy of the black hole: 



1 + - 



-|l/2~ 



(1) 



where A/bh is the black hole mass and a is the dimen- 
sionless rotation parameter (a = Jc/M^^G). 

At least two competing theories exist, however, to ex- 
plain how the system can channel this energy to pro- 
duce the collimated, relativistic jets we observe from 
GRBs. First, the energy may be extracted from the 
accretion proc ess via vv annihilation (e.g. | Woosle 



1993: MacFa dven fc ^yooslev 1999: Popham et al.l 1199 
RuflFert fc Jankal 119991 : iNaravan et al.l l2001| ). At the 
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hyper-Eddington rates expected for collapsars, the ac- 
cretion disk formed around the rotating black hole will 
be optically thick and photons are unable to escape. The 
viscous heat is instead balanced by cooling via neutrino 
emission. Annihilation of vv pairs will produce a gas of 
hot electron-positron pairs, which, assuming sufficiently 
low baryon loading, can then rapidly expand in the low- 
dens ity regions along the axis of rotation (a "fireball", 
e.g., iZhang et al.|[20Ci3l ). 

Alternatively, if the black hole is rotating and the ac- 
cretion disk is threaded by sufficiently large magnetic 
fields {B w 10^^ G for a 3M0 black hole), energy can be 
extracted directly from the rotatin g black hole via the 
Blandford-Znajek me chanism (e.g, iBlandford fc Znaiekl 
ll977HLee et al.ll2000[ ). The accretion disk serves primar- 
ily to anchor the large magnetic field (which would oth- 
erwise disperse), and the energy is emitted as a largely 
Foynting flux-dominated outflow. 

The efficiency of converting the potential energy of 
the rotating black hole plus accretion disk system into 
a form suitable to launch a collimated, relativistic out- 
flow has been intensively studied in the last decade. 
In the case of vv annihilation, the overall efficiency 
of converting the neutrino luminosity of the cooling 
disk {L^) into e+-e~ pairs (Li/p) depends sensitively on 
the mass accretion rate, the disk geometry, and (per- 
haps) the effects of general relativity on neutrino physics 
(e.g., [Pophani et al...l999; .MacFadven fc Wooslcv.1999; 



iRuffert fc Jankalll999D . But the majority of recent sim- 
ulations suggest at most a modest efficiency {rj^i, = 
LyjjjLv < 0.5%) for the energy a vailable t o launch a rel- 
ativistic jet (e g.. iRosswog et al.1 '2003: Bi rkl et al.l [20071 : 
iHarikae et al.l l2010t ). For typical disk luminosities of 
Ly K, 10^^-10^^ ergs~^, this would require an extremely 
long-lived phase of continuous accretion to explain hyper- 
energetic GRBs. 

Jets powered by MHD processes, on the other hand, 
can more easily accommodate large energy releases. To 
begin with, by extracting energy directly from the ro- 
tation of the black hole, the Blandford-Znajek pro- 
cess begins with a significantly larger energy reservoir: 
E\o\. ~ 8 X 10^"^ erg for a rapidly spinning (a — 0.9) 
black hole of mass 3 M© (Eqn. [ij . Both analytic mod- 
els (e.g.,|Lee ct al. 2000) and numerical simulations (e.g., 
lMcKinnevll2005ll suggest that as much as 5-10% of this 
rotational energy can be made available to power a colli- 
mated outflow, easily within the requirements of hyper- 
energetic events. We caution, however, that the jet ef- 
flciency is a strong function of the black hole spin, and 
slowly spinning black holes have significantly lower effi- 
ciencies. 

Setting energetics considerations aside for the moment, 
there still remains the question of whether either the vv 
annihilation mechanism or MHD processes are capable 
of producing jets with the Lorentz factors and degree of 
collimation inferred from GRB observations. Here our 
sample of i^ermi-LAT events, with their extreme initial 
Lorentz factors, offers a distinct advantage over previous 
GRB studies. Combining the Lorentz-factor limits for 
the most relativistic GRBs with the inferred jet opening 
angles from our broadband afterglow models can provide 
critical diagnostics of the jet acceleration mechanism. 

The results from MHD simulations of jet acceleration 
appear to depend sensitively on the nature of the medium 
into w hich the jet propagates. iTchekhovskov et all 
(|2009bH have recently conducted fully relativistic MHD 
simulations in which a mildly magnetized (magnetiza- 
tion parameter cr ^ 1) jet is initially confined by the 
pressure of a stellar envelope out to some radius, and 
the n allowed to propagate freely in all directions (see 
also iKomissarov et al.l[2009l ). Both authors find that the 
jet accelerates rapidly in this transition region ("rarefac- 
tion" acceleration) and can reach Lorentz factors Fq of a 
few hundred, yet still remains highly collimated (^ < 5°) 
even after leaving the region of confinement. Further- 
more, the deconfinement radii required to produce these 
outflows agree well with the expected value for Wolf- 
Rayet stars (r w 10M0"cm). 

Based on both results from these simulations and ana- 
lytical arguments, Tchekhovskoy ct al. (2009b) demon- 
strate a relationship between the initial Lorentz fac- 
tor and jet opening angle: Y^Q k. 10-30. This re- 
sult is able to nicely reproduce the typical observed 
Lore ntz fac tors (Lit hwick fc Sarill2001h and opening an- 
gles (jZeh e t al. 2003) of previous GRB samples. How- 
ever, for the events studied here, both GRB090902B 
{T^d « 70) and GRB090926A {T^d « 90) appear in- 
consistent with this result. 

One way to circumvent this requirement is if the jets 
ha ve extremely high magneti zation parameters (a ^ 
1; ITchekhovskov et all l2009al ). as such Foynting flux- 
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Fig. 10. — Two-dimensional relativistic energy release {E^^i Ej + Ske) from GRBs. Cosmologically distant {z > 0.5) events from 
the pie-Swift era are shown in red. The logarithmic mean for these events, (-Erol) = 2 X 10^^ erg, is indicated by the solid black line. 
Shaded regions correspond to Icr, 2a, and 3cr errors on this mean value. The three most nearby events (GRBs 980425, 031203, and 
060218) are plotted in green and are underluminous by several orders of magnitude. The four LAT events from this work are plotted 
in blue; all but GRB 090328 fall at the high end of the pre-Swift distribution (note that we have not plotted horizontal errors bars for 
GRB090926A due to the large uncertainty in Bke)- Instead, they are more consistent with some of the brightest events from the Swift 
era (black squares). The total relativistic energy release from GRB090926A appears to exceed 10^^ erg. Such hyper-energetic events pose 
a s evere challenge to the magnetar models, where the total energy release cannot exceed 3 X 10^^ erg (dashed black line). References — 
[Panaitescu & Kumar (2002): GRBs 990123 , 990510, 9 91208, 991216, 000301C, 010222; Yost et al. (2003): GRBs 970508, 980703, 000926; 
[Berger et al. (2004): GRBs 970508, 980703; Chevalier et al. (2004): GRB 020405: iBerger et al.l (120011): GRB 000418: iLi fc Chevalieri (il999l'l: 
GRB 980425: Soderbcre et al. (2004): GRB 031203; Soderberg et al. (2006): GRB 060218: IFrail et al.l 1(2009) : GRB 050904: IChandra et al.l 
(12008) : GRB 070125: .Cenko et al.i (,201a') : GRBs050820A, 060418, 080319B. 



dominated jets can be accelerated to extreme Lorentz 
factors over relatively large angles. However, it is unclear 
how such outflows can convert sufficient electromagnetic 
energy to accelerate electrons and produce the observed 
prompt gamma-ray emission. Other particle-acceleration 
mechanis ms besides MHD sh ocks may be required in this 
case fe. g. ■ lBeloborodovl2009( ) . Alt ernatively, the gamma - 
ray emission may be patchy (e.g.. lKumar fc Piranl j200Cl( ) 
or the jet may be structured (see, e.g., [Granot 2003 and 
references therein), so that we are measuring only the 
extrema of Fq and not the true bulk of the relativistic 
flow carrying most of the energy. 

4.3. Comparison with Other Work 

In this section, we attempt to place this work in con- 
text, both by comparing our results with those of other 
authors who have studied these same events, and by high- 
lighting additional differences between our LAT sample 
and G RBs detect ed by satellites at lower energies. 

M cBreen et al.l (pOlol) present optical and NIR obser- 
vations of three events from this work (GRBs 090323, 
090328, and 90902B), take n prim arily with the GROND 
instrument (|Greiner et al.l I2008D . In the case of 
GRB 090323, these authors find an optical spectral in- 



dex of ao — 1-90 ± 0.01, consistent with the value de- 
rived here, and our measurements of the host-galaxy flux 
agree nicely. We derive a slightly steeper optical spectral 
index, but our res ults are consist e nt at the 2.5(T level. 
Most importantly, ^Mc Breen et al.l (j2010D infer from the 
steep optical decay that a jet break occurred before the 
first optical observations began (tj < Iday), resulting 
in a narrow beaming angle {9 < 2°) and a correspond- 
ingly small coUimation-corrected prompt energy release 
[Ej < 3(1) X 10^^ erg for a constant-density (wind-like) 
circumburst medium]. We consider this possibility un- 
likely, however, as it is difficult to explain both the flat 
radio light curve and the more slowly fading X-ray after- 
glow {ax ~ 1.5) through post jet-break evolution (see 
also §1131). 

Our results also differ from those of iMcBreen et al.l 
([2010') regarding the jet break time of GRB 090328. 
These authors argue that the steep optical decay index 
{ao ~ 2.3) derived at early times requires a jet break be- 
fore the commencement of observations {tj < 1.5 days). 
The data presented in our work are not sufficient to 
uniquely determine the optical temporal decay index 
(particularly given the possible contribution from an un- 
derlying host galaxy). However, it is again difficult to 
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reconcile the shallower X-ray decay (ax ~ 1.65) and in 
particular the rising radio light curve at this time with 
models of post-jet break evolution (see also 14.41) . 

We find very similar results to those of McBrcc n et al.l 
([2010) for the temporal and spectral indices, inferred cir- 
cumburst medium, and jet break time for GRB 090902B. 
In particular, the late-time {t « 23 days) optical observa- 
tion with the VLT provides strong confirmation that the 
radio decline at t > 6 days is due to a jet break. How- 
ever, our results diverge when translating the jet break 
time into an opening angle. Because of the low density 
we derive for the circumburst medium of GRB090902B, 
a jet break time of tj w 6 days corresponds to a relatively 
narrow opening angle. The resulting constraints on 
are roughly a factor of 4 less (note that par t of this dif- 
feren ce is also due to the wider bandpass Mc Breen et al.l 
120101 use to calculate -E-v ^ iso ) . 

IPandev et al.l (|2 010f ) also present optical and X-ray ob- 
servations of GRB 090902B, most of which were included 
in our broadband modeling here. It is not surprising, 
then, that we derive similar values for the optical tem- 
poral and spectral indices (the same also holds for the 
X-ray bandpass). The afterglow model derived by these 
authors is broadly similar to ours, but we favor a some- 
what steepe r electron index (p = 2.2 vs. 1.8). Most 
importantly, IPandev et al.l (|2010D derive a limit on the 
opening angle of 6* > 6° (based on a lower limit to the jet 
break time of > 6 days). Though our derived jet break 
time is consistent with the results of these authors, our 
inferred opening angle is a factor of two smaller due to 
the lower value we have used for the circumburst den- 
sity. As a result, the limits we derive for the prompt 
g amma-ray energy release are a factor of a few less. 

IKumar fc DuranI (|2009aD ( see also IKumar fc DuranI 
l2009bD have also presented an analysis of the broadband 
afterglow of GRB090902B. Interestingly, these authors 
suggest that the delayed high-energy {E > 100 MeV) 
component observed from several LAT events is in fact 
due to the same source as the late-time afterglow emis- 
sion: synchrotron radiation from accelerated electrons 
in the circumburst medium shock-heated by the outgo- 
ing blast wave (i.e., external shock emission; see also 
iGhisellini et"aIll20Tol) . 

The afterglow parameters derived by these authors for 
GRB 090902B differ somewhat from ours, largely due 
to the fact that they assume the X-ray bandpass at 
t ^ 1 day falls below the synchrotron cooling frequency 
(i.e., i^x < Vc). While the observed X-ray spectral and 
temporal indices at i ~ 1 day are consistent with this pic- 
ture, we find that a better fit can be achieved for both 
the X-ray and optical data ii vx > Vc i ~ ^ day (see 
also lPandev"et al. 2010. 

Using the best-fit late-time afterglow parameters we 
have derived for GRB090902B (Table (H), we have at- 
tempted to calculate the magnitude of the afterglow (i.e., 
external shock) flux sX v = 100 MeV and t = 50 s. We 
find f u » 30 nJy, roughly a factor of 7 below the observed 
value (jAbdo et al.ll2009al) . The d iscrcpancv between our 
results and those of IKumar fc Du ran (2009a) is largely 
due to the lower kinetic energy we have inferred for the 
blast wave (- EKE.iso ~ 7 x 10^"^, a fact or of 18 lower than 
that used bv IKumar fc Duranll2009a|) . This would sug- 
gest that the delayed high-energy component is not due 
to afterglow emission. We caution, however, that our 
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Fig. 11. — Circumburst density and isotropic afterglow energy re- 
lease from cosmological (z > 0.5) GRBs. Constant-density models 
are shown as squares, while wind-like models are plotted as circles. 
The GRBs detected by the Fermi-LAT (red) clearly have prefer- 
entially larger isotropic energy releases and smaller circumbu rst 
densities than the rest of the sample. See the caption of Figure [TOl 
for references. 

fiux calculation does not incorporate the reduced cross- 
section for inverse Compton emission due to the Klein- 
Nishina effect at very high energies, and this will affect 
our fiux calculation to some extent. 

A further in t eresting claim from the work of 
IKumar fc DuranI (|2009af l is that the inferred value of 
the magnetic field for GRB 090902B (along with several 
other LAT-detected events) is consistent with shock com- 
pression of a modest circumstellar field {B > 30 /iG). In 
other words, no dynamo process is necessary to gener- 
ate the magnetic field strengths needed to produce the 
observed synchrotron afterglow emission. For a constant- 
density circumburst medium, the preshock magnetic field 
is given by 

B = {2TTmp€Bn)^^^c. (2) 

We also find that the product x n (and hence the de- 
rived B field) is smaller for LAT events than for other 
GRBs detected at lower energies, primarily due to the 
lower inferred circumburst densities (see below). How- 
ever, using our best-fit parameters for GRB 090902B, we 
find B 100 mG, bro adly con sisten t with the results of 
iPiran fc Nakail ()2010|) andO ()2010D . While this value is 
smaller than the preshock magnetic fields previously in- 
ferred for GRBs {B ss 10 mG for es = 0.01, n = Icm'^), 
it suggests that at least some magnetic field amplification 
may still be required. 

Finally, we return to the issue of the relatively low 
circumburst densities derived in § [31 Given the strong 
evidence for a connection between long-duration GRBs 
and broad-lined SNe Ib/c, it would be natural to expect 
the massive-star progenitors of GRBs to explode in the 
dense molecular cloud environments where we observe 
stars forming in our own Galaxy. 

Past modeling of broadband GRB afterglows, how- 
ever, has not always revealed this to be the case. Of 
the past GRBs with sufficient radio observations to 
estimate the circumburst density, the derived values 
span a large range from 1.9 x 10~^cm^^ (GRB 990123; 
iPanaitescu fc Kumar|[200l) to 680 cm^^ (GRB 050904; 
iFraU et al...2006l) . We plot' in Figure [11] the derived cir- 
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cumburst densities for these previous events (either n or 
A*, depending on the circumburst density profile) as a 
function of isotropic afterglow kinetic energy, along with 
three events from our LAT sample (we have not included 
GRB 090926A as lack of radio coverage makes density 
estimates highly degenerate). It is clear that the events 
detected by the LAT have on average larger isotropic 
energies and smaller densities than the previous sample. 

The larger isotropic kinetic energy is not hard to un- 
derstand, as the LAT is less sensitive than, for example, 
the Burst Alert Telescope (BAT; Barthclmy et al. 2005) 
onboard Swift and should therefore select brighter events 
(in terms of the high-energy fluence). Unless E^^iso and 
-EkEjIso were anticorrelated, we would expect LAT events 
to have larger isotropic blast-wave energies as well. But 
the lower densities require an additional explanation. We 
speculate on possible reasons for this difference in § [5l 

The large ratio of EKE.iso/n (or, alternatively 
i'KE,iso/A*) we have derived for these events has two 
important implications. First, the opening angle derived 
for a given jet break time scales inversely as i^KE.iso and 
proportionally to n/A^, (§ I4.4p . Thus, for a given jet 
break time, the large ratios we have inferred for LAT 
events imply smaller opening angles than were one to 
simply use a canonical value (e.g., n = 0.1 cm~'^ from 
iFrail et al.ll200ll : n = 10 cm'^ from lBloom et al.|[200l . 
In this manner, LAT events can have relatively late jet 
breaks but still be narrowly coUimated (9 < 10°). 

Secondly, a large £'KE,iso/'T- ratio will act to delay the 
deceleration time of the outflow. Assuming the "thin 
shell" case, th e outgoing relativistic blast wave will de - 
celerate when (jSari fc Piranlll999t iMolinari et "alll2007| ) 



idee (ISM): 



tdcc(Wind): 



(1 + ^) 



32TTmpC^nTQ 



f 3-Eke 



£^KE,iso(l + z) 



1/3 



(3) 

(4) 



(Not e that we have corrected Eqn. 2 from IMolinari et"al] 
120071 to remove the erroneous factor of Wp.) To some 
extent, the large ratio of EKE.iso/n will offset the ef- 
fect of the larger initial Lorentz factors for Fermi-hAT 
events. Using the parameters we have derived from 
the LAT events in our sample, we find tdoc ~ 2 s for 
GRB 090323, idee « 3s for GRB 090328, and W ~ 33s 
for GRB090902B. If the late-time GeV emission is in- 
deed due to external shock emission, the delay between 
the MeV and GeV photons should correspond roughly 
to the blast-wave deceleration time. This may be prob- 
lematic for GRB090902B, for which the observed delay 
(~ 4 s) is significantly less than our inferred deceleration 
time. We caution, however, that the large systematic 
uncertainties in £'KE,iso, n, and particularly Fq (§ 14. 4|) 
may introduce a significant error into calculations of the 
deceleration time. 

4.4. Limitations and Future Work 

To better understand our results in § 14.11 § 14.21 and 
§ 14.31 on GRB energetics and circumburst density, it is 
important to consider the robustness of our modeling. 
The most crucial factor in determining the true energy 
release is the coUimation correction, which we therefore 
examine in more detail here. 



For a wind-like circumburst medium , the jet half- 
open ing angle 9 can be written as (e.g., iChevalier fc Lil 
[2OOOI) 



6l(Wind) =0.10 



l + z 



-1/4 



/ Eke,!'. 



A, 



1.0 cm 



\ 10^^ erg 

1/4 / ^. X 1/4 



-1/4 



^ ) [rad]. (5) 



This result assumes only that the shock is ul- 
trarelativistic and underg oes a self-similar evolution 
(|Blandford fc McKed 119761 ). and is therefore relatively 
robust. 

Only in the case of GRB090926A (and also for 
GRB 090 902B if we include th e late-time VLT observa- 
tion from lMcBreen et al.ll20Tol ) do we find clear evidence 
for a jet break in multiple bandpasses. For GRB 090323, 
the lack of a steepening in the X-ray and optical light 
curves apparently limits tj > 10 days. If we simply ig- 
nore our modeling results and assume that the shock 
energy is comparable to the prompt gamma-ray energy 
(^-KE.iso ~ 4 X 10'^'* erg) and that the progenitor wind 
speed and mass-loss rate are comparable to those ob- 
served from Galactic Wolf-Rayet stars (A^, « Igcm^^), 
then we find a lower limit on the opening angle of 
6'(Wind) > 3°. The resulting limit on the prompt en- 
ergy release is therefore i?^ > 5 x 10^^ erg, similar to 
what we derived from our broadband models. 

If we were instead to assume the GRB exploded in 
a constant- density medium, the opening angle then be- 
comes (e.g.. lSari et al.llT999( ) 



6'(ISM) =0.12 



1 



Acm~ 



\ 10^^ erg 



-1/8 



3/8 



adl 



(6) 



Under similar assumptions as above, the limit on the 
opening angle would be 0(ISM) > 6°, corresponding to 
a prompt energy release oi E^ > 2 x 10^^ erg. Clearly 
a constant-density circumburst medium only makes the 
energy requirements more strict. 

An alternative possibility for GRB 090323 
()McBreen et al.l 120101 : also suggested by iSchadv et al.l 
[20071 to explain the afterglow of GRB 061007) is a 
jet break before the beginning X-ray and optical 
observations {t < 1 day for GRB 090323). Much like 
GRB 061007, however, the implied post-break decay 
index in the X-ray {ax ~ 1-5) is too shallow to accom- 
modate a typical electro n spectral index p (for 1/ > Vm, 
fy oc t~P post iet-break: iSari et allll999| ). The observed 
behavior would therefore require late-time energy input 
to flatten the X-ray and optical decay. Furthermore, 
the relatively flat radio light curve for t <20 days would 
require that the synchrotron self-absorption frequency 
fell above the radio bandpass during this phase; oth- 
erwise the radio bandpass would be required to decay, 
either as t^^/^ [v < v^) or t^^ {v > J^m 0- We consider 

^'^ We note that the above results are independent of circumburst 
medium, as the expansion of the outflow is predominantly lateral 
after the jet break ( Chevalier fc Li 2 000 ). 
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this somewhat contrived picture unhkely in the case of 
GRB 090323, particularly given the reasonable quality 
of our fits using only standard afterglow theory. 

A similar analysis of GRB 090328 suggests that the 
jet-break time cannot occur any earlier than t « 4 days, 
based again on the relatively steady decay seen in the 
X-ray light curve and, in particular, on the rising radio 
emission at early times. The corresponding limits on 
the opening angles are therefore ^(Wind) > 8°, ^(ISM) 
> 7° . Both values result in a prompt energy release Ej > 
10^^ erg, broadly consistent with our modeling results. 

In the end, all of these estimates of energy, den- 
sity, and jet geometry rest on the relativistic syn- 
chrotron model. While the standard afterglow model 
has undergone continuous improvements and it has 
been well tested (Mcszaros & Recs 1997; Sa ri et al. 
19991: IWiiers fc Galamai Il999l: iGhevalier fc Li 
KTaitescu fc Kumad 120001 : lYost et all 120031 : 



1998, 



200C; 



Willingale et al.ll2007D . it still must make several simpli- 



fying assumptions about the shock dynamics, magnetic 
field generation, particle acceleration, energy injection, 
and the circumburst density structure. Of particular con- 
cern are the recent results of relativistic simulations by 
iZhang fc MacFadvenI (|2009| ). which suggest that viewing 
angle can dramatically alter the observed afterglow emis- 
sion at early times. 

If we are to verify this class of hyper-energetic bursts, 
we need independent est i mate s of the relativistic energy 
content. iWaxman et al.l ()1998[ ) first suggested that late- 
time (radio) calorimetry could be used for this purpose 
in order to sidestep early-time complications such as the 
outflow geometry, the density structure, and ongoing ac- 
tivity from the central engine. Recent relativistic hy- 
drodynamic simulations show that the earlier analytic 
models which used spherical geometry and Sedov- Taylor 
dynamics to explain the late-time behavior were approxi- 
mately correct and produce a robust estima te of the total 
kinetic energy ()Zhang fc MacFadven|[2009l ). In the past, 
the method has been limi ted to only a small numb er of 
radio-bright events (e.g., Ivan der Horst et al.l 120081 and 
references therein). But it is somewhat encouraging that 
both early- and late-time estimates for the most ener- 
getic pre- Fermi eve nt (GRB 970508; ^ke » lO^'^ erg) ar e 
in close agreement (|Yost et al.ll2003t iBerger et al.ll2004[ ). 
With the advent of radio facilities having greatly in- 
creased sensitivi ty such as LOFAR (|van Haarleml[2nn5h 
and the EVLA (|Perlev et al.ll2009d) . it should soon be 
possible to verify candidate hyper-energetic events by 
searching for long-lived {t > 100 days) radio afterglows. 

Finally, we add a further note of caution that the 
derivation of the initial Lorentz factor Fq is also 
subject to significa nt systematic uncertainties (e.g., 
iBosniak et al.ll2009t ). In order to calculate the optical 
depth to pair production, it is necessary to determine 
the effective blast-wave radius as well as the instanta- 
neous spectral parameters. The former is typically in- 
ferred from the variability time scale, while the latter is 
usually averaged over the entire duration of the prompt 
emission, and it is not always clear how large an effect 
this will have on the calculation of Fq. The calculations 
are further complicated if the spectrum is not a pure 
power law, as would be the ca se for thermal emi ssion aris- 
ing from the jet photosphere ()Rvde et al.ll2010|) . Here we 



have attempted to calculate the initial Lorentz factors in 
a consistent manner, and it is not unreasonable to believe 
that the relative ordering of Fq is robust. However, small 
changes in Fq can have a large effect on some of the pa- 
rameters derived here (e.g., ^dcc; Eqns.[3]and|4l). Future 
early afterglow observations of Fermi-LAT GRBs to di- 
rectly m easure the deceleratio n time and hence constrain 
Fq (e.g.. iMolinari et al]|2007[ ) could provide valuable in- 
sight in this area. 

5. CONCLUSIONS 

We have undertaken extensive broadband continuum 
(radio, optical, and X-ray) and spectroscopic observa- 
tions of four long-duration GRBs (GRBs 090323, 090328, 
090902B, and 090926 A) detected by the LAT instru- 
ment on the Fermi satellite at GeV energies. This work 
was motivated by the realization that Fermi is espe- 
cially sensitive to GRBs with large isotropic energy re- 
lease, and hence provides an interesting sample of events 
to test GRB central-engine models and their relativis- 
tic outflows. Our afterglow models constrain the jet 
break times and the density of the circumburst medium, 
from which we derive the collimation angle and hence 
beaming-corrected energy release for each event. We find 
three GRBs with a total relativistic content about an 
order of magnitude in excess of the canonical 10^^ erg, 
with GRB 090926A almost certainly in excess of 10^^ erg. 
This analysis provides support for our earlier claim of a 
class of hyper-energetic GRBs. The discovery of more 
GRBs with total energy release > 10^^ erg is troubling 
for central engines in which the energy to drive the jet 
is derived either from a rotating magnetar or collapsars 
powered by neutrino annihilation. For this reason, we are 
led to believe that, at least for hyper-energetic GRBs, the 
massive star progenitor collapses directly to a black hole 
and the rotational energy of this system is extracted via 
the Blandford-Znajek process. 

Although we find relatively narrow opening angles for 
all four events {9 < 10°), the extreme initial Lorentz fac- 
tors inferred for these LAT events imply that the product 
OTq can be a factor of 5-10 larger than estimates of previ- 
ous GRBs detected at MeV energies. These values are in- 
consistent with recent simulations of low-magnetization 
MHD jets, suggesting that the outflow may be at least 
initially Poynting-fiux dominated. If this is indeed the 
case, it is unclear how the initial kinetic energy of the 
outflow is converted to prompt gamma-ray emission. 

Interestingly, for the three events having sufficient ra- 
dio coverage to derive a circumburst density, we find 
anomalously large values of EKE.iso/n (or, for a wind- 
like medium, i?KE,iso/^*)- While the large £'KE,iso values 
are simple to understand, the low circumburst densities 
require a more complex explanation. 

One possibility is that the progenitor stars of LAT 
GRBs are somehow different from the progenitors of most 
previous GRBs detected at MeV or keV energies. It 
is currently thought that GRB progenitors are distin- 
guished from the progenitors of ord inary SN e Ib/c by 
their low metallici ties (e . g.. MacFadven fc ^Yoos lcv 199a; 
iWooslev fc ifegeH l200a iModiaz et al.ll2008f ): the lower 
mass-loss rates allow the progenitors of GRBs to keep 
more angular momentum. The increased rotation evac- 
uates a cavity through which a relativistic jet can prop- 
agate. If LAT events have larger initial Lorentz factors. 
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it may be that they come from lower metalhcity progeni- 
tors with minimal pre-explosion mass loss. Observations 
of the host galaxies of these events, both through absorp- 
tion and emission spectroscopy, may help shed light on 
this matter. 

It is also possible that the low density preference is the 
result of other, more subtle, selection effects. In particu- 
lar, if the GeV emission arises in the external shock, the 
Fermi-ljAT could be biased towards events in low-density 
environments. If the circumburst density is too high, the 
blastwave will decelerate at small radii (Equations [3] and 
0]), where the outflow may be opaque to GeV photons. 
More observations of LAT GRBs, particularly at very 
early times, would help to investigate this hypothesis. 

We end by emphasizing the importance of afterglow 
observations of high-£^^_iso events in the Fermi era to 
provide further confirmation of this picture. Such GRBs 
are either highly collimated outflows {9 <2°) with a typ- 
ical energy release, or truly hyper-energetic events; both 
represent extreme tests of jet collimation and central- 
engine models, respectively. Current efforts suffer from 
delays in LAT localizations and limited ground-based af- 
terglow follow-up efforts. The latter can be improved by 
fo cusing rare follow-up reso urces on Fermi-L AT GRBs ; 
as iNvsewander et al.l (|2009D and iMcBreen et~all (|20100 
have shown, these events have brighter X-ray and op- 
tical afterglows on average, and are therefore accessible 
even for moderate-aperture optical facilities. Targeting 
these bright afterglows will make it easier to measure the 
jet breaks, which have proven almost impossible to ob- 
tain in the Swift era. Finally, we note that one testable 
consequence of hyper-energetic GRBs is long-lived after- 
glow emission (> lyr). If the shock microphysics and 
the circumburst density do not undergo a drastic evolu- 
tion, it should be possible to detect these afterglows with 
the upcoming generation of radio facilities and carry out 



calorimetry measurements. 
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^ UT at beginning of exposure. 

^ Time from midpoint of exposure to i^ermz-GBM trigger. 

^ Reported magnitudes have not been corrected for Galactic extinction {E(B — V) — 0.025 mag JSchlegel et aLlUQQ*^ . Obser- 
vations in the R band arc referenced to Vega, while all other filters arc reported on the AB magnitude system jOke & GunnI 
[1981). 

^ * - This work; 1 -lUpdikc ct al. (2009b); 2 - Wang ct al. (2009); 3 - Kann et al. (20093); 4 -[Burenin ct'aL ("20W); 5 - 
Perlev et al.| J2009a|); 6 - Guidorzi ct al. (2009); 7 - Pcrlcy (2009); 8 - Kann ct al.. (.2009*3) : 9 - [Kann ct al. (.2009a.); 10 - 
Rumvantsev & Pozanenkd fl2009l '); 11 - |de~Ugarte Postigo e^t"a,n l l2009j) . 
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TABLE 10 

Optical/NIR Observations of GRB 090328 
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19.54 ±0.03 


1 


2009 Mar 29.98 


1.60 


GROND 


J 


480.0 


19.54 ±0.06 


1 


2009 Mar 29.98 


1.60 


GROND 


H 


480.0 


19.02 ±0.06 


1 


2009 Mar 29.98 


1.60 


GROND 


K 


480.0 


18.52 ±0.08 


1 


2009 Mar 30.06 


1.69 


UVOT 


U 


4111.2 


20.80 ±0.11 




2009 Mar 31.42 


3.04 


UVOT 


U 


4162.2 


21.96 ±0.20 




2009 Apr 3.79 


6.40 


UVOT 


u 


2094.3 


21.52 ±0.24 




2009 Apr 4.79 


7.40 


UVOT 


u 


2507.4 


21.69 ±0.23 




2009 Apr 9.68 


12.33 


UVOT 


u 


9493.6 


22.67 ±0.24 





^ UT at beginning of exposure. 

^ Time from midpoint of exposure to Fermi-G^M trigger. 

Reported magnitudes have not been corrected for Galactic extinction {E(B — V) — 0.057 mag MSchlegel et aLlUQQ^ . Obser- 
vations in the U band are referenced to Vega, while all other filters are reported on the AB magnitude system l lOke & GunnI 
[1981). 

* — This work; 1 - lUpdike eTaT] j2009al '). 



TABLE 11 
Optical Observations of GRB 090902B 



Date'' 
(UT) 


Time Since Burst'' 
(days) 


Telescope/Instrument 


Filter 


Exposure Time 
(s) 


Magnitude'^ 


Reference'* 


2009 Sep 2.51 


0.06 


ROTSE-IIIa 


R 


890.0 


16.4 ±0.5 


1 


2009 Sep 2.72 


0.27 


ROTSE-IIId 


R 


846.0 


> 18.7 


1 


2009 Sep 2.98 


0.53 


UVOT 


U 


1074.7 


20.31 ±0.15 


* 


2009 Sep 3.08 


0.64 


UVOT 


U 


3612.9 


20.54 ± 0.09 


* 


2009 Sep 3.19 


0.74 


Nickel 


R 


3300.0 


20.60 ±0.10 


1 


2009 Sep 3.24 


0.79 


Gemini-N/GMOS 


r' 


180.0 


20.68 ±0.15 


* 


2009 Sep 3.24 


0.82 


UVOT 


U 


8756.0 


20.92 ±0.08 


* 


2009 Sep 3.33 


0.88 


Liverpool 


R 


1800.0 


21.04 ±0.11 


1 


2009 Sep 3.74 


1.30 


UVOT 


U 


3852.8 


22.06 ±0.26 


* 


2009 Sep 3.87 


1.42 


WHT-LIRIS 


J 


1800.0 


19.99 ±0.15 


1 


2009 Sep 3.89 


1.44 


Liverpool 


R 


1800.0 


21.40 ±0.10 


1 


2009 Sep 3.92 


1.47 


WHT-LIRIS 


K 


2592.0 


18.92 ±0.20 


1 


2009 Sep 3.99 


1.54 


GROND 


r' 


738.0 


21.54 ±0.05 


1 


2009 Sep 4.00 


1.55 


NOT 


V 


900.0 


21.67 ±0.11 


1 


2009 Sep 4.01 


1.56 


NOT 


R 


900.0 


21.40 ±0.11 


1 


2009 Sep 4.02 


1.57 


NOT 


I 


900.0 


20.72 ±0.11 


1 


2009 Sep 4.36 


1.91 


UKIRT-WFCAM 


J 


1080.0 


20.20 ±0.20 


1 


2009 Sep 4.36 


1.91 


UKIRT-WFCAM 


K 


1080.0 


18.90 ±0.25 


1 


2009 Sep 4.98 


2.53 


GROND 


r' 


738.0 


22.01 ±0.07 


1 


2009 Sep 5.93 


3.48 


Liverpool 


R 


1800.0 


22.60 ±0.25 


1 


2009 Sep 6.37 


3.97 


UVOT 


U 


10549.1 


22.25 ±0.21 


* 


2009 Sep 8.98 


6.53 


GROND 


r' 


738.0 


23.07 ±0.17 


1 


2009 Sep 10.37 


7.97 


UVOT 


U 


9969.0 


> 22.56 


* 



^ UT at beginning of exposure. 

^ Time from midpoint of exposure to i^ermi-GBM trigger. 

^ Reported magnitudes have not been corrected for Galactic extinction {E{B — V) — 0.042 mag; ISchlegel et al.||199SD . 
Observations in the the r' filter are reported on the AB magnitude system ]Oke Gunn|[l983l) , while all other filters are 
referenced to Vega. 

<i * - This work; 1 - iPandev eTaP EoTdi ). 



Afterglow Observations of FermirLAT GRBs 



TABLE 12 

Optical/NIR observations of GRB090926A 



Date'' 
(UT) 


lime Since Burst 


Telescope/Instrument 


Filter 
(s) 


Exposure Time 


Magnitude^ 


2009 Sep 26.72 


0.5415 


UVOT 


u 


358.0 


18 33 -1- 06 


200Q Son 2fi 72 


0.5501 


UVOT 


V 


358.0 


XO.Orl _1_ VJ. XW 


200Q Srn 2fi 78 


5989 


UVOT 


If 


675.0 


xo.iJU _i_ yj.yjtj 


2009 Son 2fi 79 


0.6149 


UVOT 


V 


675.0 


1201 -1- n no 

xo.c/x _i_ yj.yjc/ 


2009 Ron 2fi 84 


0.6649 


UVOT 


If 


704.5 


1 2 71 -L QK 
xo. 1 X _i_ yj.yjtj 


2009 Spd 2fi 8^ 


0.6816 


UVOT 


V 


703.3 


1 Q no -1- n 1 n 

XCJ.VC/ _1_ \J.±\J 


2009 Sen 2fi 90 


0.7299 


UVOT 


u 


526.0 


18 84 -1- 07 

X . yJ^ _l \Jm\J i 


2009 Sen 2fi 9fi 


0.7966 


PROMPT-2 


V 


560.0 


1 a 7R + n 22 


2009 Sen 2fi 9fi 


0.7969 


PROMPT-4 


R 


640.0 


1 8 69 -1- 1 6 

XO.Ui? _1_ L/.X(J 


2009 Sen 2fi 9fi 


0.7986 


PROMPT-5 


J 


400.0 


1 Q no -1- n 1 Q 

XO.Uii _1_ \J.XC7 


2009 Sen 26 98 


8062 


UVOT 


u 


459.4 


1 8 94 -1- 08 


2009 Son 26 99 


0.8172 


UVOT 


V 


459.4 


18 97 -1- 1 2 


2009 Sep 27.00 


0.8250 


PROMPT-5 


I 


720.0 


18.12 ± 0.07 


2009 Sen 27 00 


8259 


PROMPT-4 


R 


880.0 


1 8 49 -1- 05 


2009 Sen 27 00 


0.8314 


PROMPT-2 


V 


1600.0 


19 12-1-0 10 


2009 Sep 27.00 


0.8372 


PROMPT-5 


I 


880.0 


18.06 ± 0.05 


2009 Sep 27.00 


0.8379 


PROMPT-4 


R 


880.0 


18.45 ± 0.04 


2009 Sep 27.00 


0.8553 


PROMPT-4 


R 


800.0 


18.41 ± 0.04 


2009 Sep 27.00 


0.8553 


PROMPT-5 


I 


800.0 


18.04 ± 0.05 


2009 Sep 27.00 


0.8583 


SMARTS / ANDICAM 


I 


2160.0 


18.17 ± 0.07 


2009 Sen 27 00 


8583 


SMARTS / ANDICAM 


J 


1800.0 


17 28 -1- 07 


2009 Sen 27 00 


8603 


PROMPT-2 


V 


1520.0 


19 12 -1- 09 


2009 Sen 27 04 


0.8654 


PROMPT-3 




880.0 


Xi7.(JU _1_ (J.XL» 


2009 Sen 27 04 


8659 

U.OUUi? 


PROMPT-4 


R 


800.0 


1 a QR -1- n 04 

XO.tlU _1_ (J.U^ 


2009 Sen 27 04 


8660 


PROMPT-5 


J 


800.0 


1 7 97 -1- 05 

X 1 . 1 _1_ \J . \JyJ 


2009 Sen 27 04 


8662 


UVOT 


If 


689.8 


1 2 an -1- n 06 

XO.OU _1_ u.uu 


2009 Sen 27 04 


0.8761 


PROMPT-4 


R 


800.0 


1 2 QR -1- n 04 

XO.OU _1_ \J.yj'-t 


2009 Sen 27 04 


0.8764 


PROMPT-5 


J 


640.0 


1 7 qc; -1- n 05 


2009 Sen 27 04 


8808 

U . OO vy O 


PROMPT-2 


V 


1520.0 


19 19 -1- 09 


2009 Sen 27 04 


0.8811 


PROMPT-3 


B 


1520.0 


1 9 57 -1- 07 

X c • c _i_ yj *yf 1 


2009 Sen 27 0^ 


8826 

U.OOi;U 


UVOT 


V 


689.7 


1 2 SR -1- n OQ 

XO.OU _1_ yJ-yj^y 


2009 Sen 27 04 


8858 


PROMPT-5 


J 


800.0 


1 7 87 -1- 04 


2009 Sen 27 04 


0.8864 


PROMPT-4 


R 


800.0 


18 31 + 04 


2009 Sen 27 04 


8962 


PROMPT-5 


J 


800.0 


1 7 92 -1- 05 


2009 Sen 27 04 


0.8967 


PROMPT-4 


R 


800.0 


18 29 -1- 04 


2009 Sen 27 04 


0.9012 


PROMPT-3 


B 


1600.0 


19 56 -1- 07 


2009 Sen 27 04 


0.9013 


PROMPT-2 


V 


1520.0 


19 02 -1- 08 


2009 Sen 27 08 


9055 


PROMPT-5 


J 


640.0 


1 7 85 -1- 05 


2009 Sen 27 08 


9065 


PROMPT-4 


R 


720.0 


1 2 54 -1- n 04 


2009 Sen 27 08 


0.9172 


PROMPT-4 


R 


880.0 


1 2 oQ -1- n 04 

XO.iii/ _1_ yj.yj'-t 


9nnQ ^pn 97 n^^ 


0.9172 


PROMPT 5 


J 


880.0 


1 7 70 _|_ n n'i 


2009 Son 27 08 


0.9223 


PROMPT-3 


B 


1680.0 


1 Q Qo _|_ n (-15 
xc.^jo _i_ yj*yj^ 


2009 Sen 27 08 


9225 


PROMPT-2 


V 


1600.0 


1 2 qc; -1- n 07 

XO.t/tJ _1_ yJ-yJI 


2009 Sen 27 08 


9280 


PROMPT-4 


R 


800.0 


1 8 29 -1- 04 


2009 Sen 27 08 


0.9284 


PROMPT-5 


J 


560.0 


17 85-1-0 05 


2009 Sen 27 10 


0.9331 


UVOT 


u 


691.7 


1 8 64 -1- 05 


2009 Sen 27 08 


9388 


PROMPT-4 

X X \/V_y IVX X J_ rr 


R 


880.0 


1 2 1 K _i_ n 03 

XO.XfJ _1_ u.uu 


9nnQ CJfsr* 97 OS 


u.yooo 


PROMPT-5 


I 


880.0 


1 7 7fi _|_ n fiA 
X 1 . 1 u ^ u.u^ 


9nnQ Spn 97 08 


0.9419 


PROMPT-3 

X X lyV-/±VXX X_ t_l 


B 


1360.0 


1 Q cc _(_ n no 

Xi7.UU _1_ U.UO 


2009 Sen 27 08 


0.9421 


PROMPT-2 

X X J.VX X X- ^ 


V 


1280.0 


1 a at; _i_ n OS 

XO.OU _1_ U.UO 


2009 Sen 27 1 S 


0.9475 


PROMPT-4 

X X LpV^ IVX X X rr 


R 


480.0 


1 a -1- 05 

XO.iiU _1_ u.uu 


2009 Sen 27 1 S 


0.9477 


PROMPT-5 

X X IVX X A- yj 


J 


480.0 


1 7 Qo -1- n 05 

X ( .Oil _i_ u.uu 


2009 Sen 27 12 


0.9495 


UVOT 


V 


691.7 


1 2 7R -1- n 09 

XO. 1 U _1_ U.Ui? 


2009 Sen 27 1 7 


9933 


PROMPT-3 


B 


960.0 


19 60 -1- 07 

i-Cf *yjyj _i_ yj *yj \ 


2009 Sen 27 1 7 


9933 


PROMPT-5 


J 


960.0 


1 7 82 -1- 04 


2009 Sen 27 17 


9935 


PROMPT-2 


Y 


960.0 


18 88 -1- 10 

xo.oo _i_ yj * i-yj 


2009 Sen 27 17 


9938 


PROMPT-4 

X XV.V-'IVXX X ± 


R 


1040.0 


1201 -1- n 03 

xo.il X _i_ yj.yjij 


2009 Sen 27 17 


9992 


UVOT 


u 


554.0 


1 8 65 -1- 06 

xo.u^w^ _1_ yj *yjyj 


2009 Sen 27 1 7 


1.0058 


PROMPT-5 


J 


960.0 


1 7 83 -1- 03 

X i .ov* _i_ yj *yj^ 


2009 Sen 27 17 


1.0062 


PROMPT-4 

X X l/Xw* IVX X X rr 


R 


880.0 


1 2 24 -1- n 03 

XO.iirt _1_ U.UU 


2009 Sen 27 17 


1.0100 


PROMPT-3 


B 


1200.0 


19 43 -1- 06 


2009 Sen 27 17 


1.0107 


PROMPT-2 


V 


1520.0 


19 03 -1- 09 


2009 Sen 27 1 7 


1.0166 


PROMPT-4 


R 


560.0 


1 a 1 Q -1- n 04 

X . X £7 _1_ \J . \J'-t 


2009 Sep 27.17 


1.0173 


PROMPT-5 


I 


640.0 


17.79 ± 0.04 


2009 Sep 27.21 


1.0293 


PROMPT-4 


R 


960.0 


18.26 ±0.03 


2009 Sep 27.21 


1.0295 


PROMPT-5 


I 


880.0 


17.85 ±0.04 


2009 Sep 27.21 


1.0352 


PROMPT-2 


V 


1680.0 


18.89 ±0.07 


2009 Sep 27.21 


1.0363 


PROMPT-3 


B 


1360.0 


19.53 ±0.07 


2009 Sep 27.21 


1.0408 


PROMPT-5 


I 


800.0 


17.87 ±0.04 


2009 Sep 27.21 


1.0412 


PROMPT-4 


R 


880.0 


18.43 ±0.04 


2009 Sep 27.21 


1.0530 


PROMPT-5 


I 


960.0 


18.01 ± 0.04 


2009 Sep 27.21 


1.0530 


PROMPT-4 


R 


960.0 


18.42 ±0.04 


2009 Sep 27.21 


1.0590 


PROMPT-2 


V 


1840.0 


19.13 ±0.10 


2009 Sep 27.21 


1.0592 


PROMPT-3 


B 


1920.0 


19.64 ± 0.06 



Cenko et al. 



TABLE 12 — Continued 



Date^ 




J. tJltJoCO^t;/ lllDtl LillifcJIit 


Filter 




rrni+nrlo*^ 
iVlclglll t LI Lit; 




\°-) 




\°) 






2009 Sep 27.21 


1.0653 


PROMPT-4 


R 


960.0 


18.43 ± 0.04 


2009 Sep 27.21 


1.0655 


PROMPT-5 


I 


960.0 


18.00 dz 0.05 


2009 Sep 27.24 


1.0658 


UVOT 


U 


730.1 


18.67 ± 0.05 


2009 Sep 27.25 


1.0777 


PROMPT-4 


R 


960.0 


18..56 ± 0.04 


2009 Sep 27.25 


1.0780 


PROMPT-5 


I 


960.0 


18.02 ±0.04 


2009 Sep 27.25 


1.0839 


PROMPT-3 


B 


1920.0 


19.59 ± 0.05 


2009 Sep 27.25 


1.0839 


PROMPT-2 


V 


1840.0 


19.15 ±0.09 


2009 Sep 27.25 


1.0900 


PROMPT-4 


R 


960.0 


18.45 ± 0.03 


2009 Sep 27.25 


1.0902 


PROMPT-5 


I 


880.0 


17.97 ± 0.04 


2009 Sep 27.25 


1.1028 


PROMPT-5 


I 


1040.0 


17.98 ±0.04 


2009 Sep 27.25 


1.1029 


PROMPT-4 


R 


1040.0 


18.39 ± 0.03 


2009 Sep 27.25 


1.1091 


PROMPT-3 


B 


2000.0 


19.51 ±0.04 


2009 Sep 27.25 


1.1093 


PROMPT-2 


V 


1920.0 


18.96 ±0.07 


2009 Sep 27.29 


1.1158 


PROMPT-4 


R 


960.0 


18.47 ±0.03 


2009 Sep 27.29 


1.1158 


PROMPT-5 


I 


960.0 


18.04 ±0.04 


2009 Sep 27.29 


1.1283 


PROMPT-5 


I 


960.0 


18.02 ±0.04 


2009 Sep 27.29 


1.1287 


PROMPT-4 


R 


1040.0 


18.45 ±0.03 


2009 Sep 27.29 


1.1348 


PROMPT-3 


B 


2000.0 


19.62 ±0.05 


2009 Sep 27.29 


1.1352 


PROMPT-2 


V 


1920.0 


19.09 ± 0.09 


2009 Sep 27.29 


1.1413 


PROMPT-5 


I 


1040.0 


17.97 ± 0.04 


2009 Sep 27.29 


1.1415 


PROMPT-4 


R 


960.0 


18.43 ± 0.03 


2009 Sep 27.33 


1.1545 


PROMPT-4 


R 


960.0 


18.47 ±0.03 


2009 Sep 27.33 


1.1550 


PROMPT-5 


I 


1040.0 


18.06 ±0.05 


2009 Sep 27.33 


1.1617 


PROMPT-2 


V 


2000.0 


19.11 ±0.07 


2009 Sep 27.33 


1.1617 


PROMPT-3 


B 


2080.0 


19.67 ±0.06 


2009 Sep 27.33 


1.1684 


PROMPT-4 


R 


1040.0 


18.51 ±0.03 


2009 Sep 27.33 


1.1685 


PROMPT-5 


I 


1040.0 


18.06 ±0.05 


2009 Sep 27.33 


1.1818 


PROMPT-4 


R 


1040.0 


18.52 ±0.04 


2009 Sep 27.33 


1.1826 


PROMPT-5 


I 


960.0 


18.24 ±0.06 


2009 Sep 27.33 


1.1850 


PROMPT-3 


B 


1520.0 


19.74 ±0.08 


2009 Sep 27.33 


1.1859 


PROMPT-2 


V 


1600.0 


19.05 ±0.10 


2009 Sep 27.33 


1.1919 


PROMPT-4 


R 


480.0 


18.63 ± 0.07 


2009 Sep 27.33 


1.1923 


PROMPT-5 


I 


480.0 


18.28 ±0.11 


2009 Sep 27.38 


1.2093 


PROMPT-3 


B 


480.0 


19.73 ±0.19 


2009 Sep 27.38 


1.2093 


PROMPT-4 


R 


480.0 


18.71 ± 0.08 


2009 Sep 27.38 


1.2093 


PROMPT-5 


I 


480.0 


18.18 ±0.10 


2009 Sep 27.38 


1.2095 


PROMPT-2 


V 


480.0 


18.97 ± 0.23 


2009 Sep 27.84 


1.6702 


UVOT 


U 


663.2 


19.67 ±0.10 


2009 Sep 27.91 


1.7349 


UVOT 


U 


735.0 


19.74 ±0.09 


2009 Sep 27.95 


1.7846 


UVOT 


V 


2168.5 


19.82 ±0.09 


2009 Sep 27.97 


1.8018 


UVOT 


u 


737.0 


19.71 ±0.09 


2009 Sep 28.04 


1.8685 


UVOT 


u 


742.9 


19.78 ±0.10 


2009 Sep 28.11 


1.9344 


UVOT 


u 


577.3 


19.79 ±0.11 


2009 Sep 28.12 


1.9476 


UVOT 


u 


341.3 


19.73 ±0.14 


2009 Sep 28.40 


2.2343 


UVOT 


u 


1329.2 


20.13 ±0.08 


2009 Sep 28.55 


2.3829 


UVOT 


u 


1730.2 


20.01 ± 0.07 


2009 Sep 28.56 


2.3950 


UVOT 


V 


1131.6 


20.07 ±0.14 


2009 Sep 28.88 


2.7199 


UVOT 


u 


2642.7 


20.11 ±0.06 


2009 Sep 28.89 


2.7277 


UVOT 


V 


2417.0 


20.11 ±0.10 


2009 Sep 29.08 


2.9301 


PROMPT-5 


I 


3200.0 


19.48 ±0.10 


2009 Sep 29.08 


2.9346 


PROMPT-4 


R 


2560.0 


19.96 ±0.12 


2009 Sep 29.13 


2.9642 


PROMPT-5 


I 


1840.0 


19.39 ±0.11 


2009 Sep 29.13 


2.9644 


PROMPT-4 


R 


1840.0 


19.88 ±0.12 


2009 Sep 29.16 


2.9892 


SMARTS / ANDICAM 


I 


2160.0 


19.53 ±0.07 


2009 Sep 29.16 


2.9892 


SMARTS / ANDICAM 


J 


1800.0 


18.35 ±0.09 


2009 Sep 29.17 


3.0041 


PROMPT-5 


I 


2720.0 


19.47 ±0.10 


2009 Sep 29.17 


3.0055 


PROMPT-3 


B 


3520.0 


> 21.85 


2009 Sep 29.17 


3.0056 


PROMPT-4 


R 


2960.0 


19.85 ± 0.09 


2009 Sep 29.17 


3.0083 


PROMPT-2 


V 


3600.0 


> 20.94 


2009 Sep 29.53 


3.3664 


UVOT 


V 


2533.9 


20.75 ±0.13 


2009 Sep 29.54 


3.3796 


UVOT 


u 


3233.8 


20.58 ± 0.07 


2009 Sep 30.00 


3.8439 


PROMPT-5 


I 


5840.0 


19.78 ±0.11 


2009 Sep 30.00 


3.8446 


PROMPT-4 


R 


6560.0 


20.28 ±0.10 


2009 Sep 30.08 


3.9295 


PROMPT-5 


I 


6000.0 


19.84 ±0.12 


2009 Sep 30.08 


3.9302 


PROMPT-4 


R 


5680.0 


20.54 ±0.14 


2009 Sep 30.17 


4.0116 


PROMPT-3 


B 


30640.0 


21.74 ±0.13 


2009 Sep 30.17 


4.0121 


PROMPT-2 


V 


28720.0 


20.98 ±0.17 


2009 Sep 30.17 


4.0249 


PROMPT-5 


I 


7440.0 


19.99 ±0.12 


2009 Sep 30.17 


4.0253 


PROMPT-4 


R 


7520.0 


20.44 ± 0.09 


2009 Sep 30.29 


4.1428 


PROMPT-5 


I 


10400.0 


20.07 ±0.11 


2009 Sep 30.29 


4.1445 


PROMPT-4 


R 


10400.0 


20.62 ± 0.10 


2009 Sep 30.55 


4.3878 


UVOT 


V 


3091.7 


20.60 ±0.11 


2009 Sep 30.55 


4.3944 


UVOT 


U 


3341.7 


20.90 ± 0.08 


2009 Oct 1.04 


4.8842 


PROMPT-4 


R 


12640.0 


20.62 ±0.10 


2009 Oct 1.04 


4.8854 


PROMPT-5 


I 


12400.0 


20.06 ±0.10 



Afterglow Observations of Fermi-LAT GRBs 

TABLE 12 — Continued 



29 



Datc^" 


Time Since Burst'' 


Telescope /Instrument 


Filter 


Exposure Time 


Magnitude*^ 


(UT) 


(d) 




(s) 






zuuy Kjcx i.ii 


4.yDoo 


ClAAARTQ / A NTFlTr^ A 1\/T 


T 
1 


01 f^c\ n 

zlDU.U 


on _i_ n 1 n 
zu.oo it U. lU 


onno r^^.+ i 1 1 
zUUy Uct l.il 


A Ofic; Q 
4.yDOo 


oiVlArvi / AlNJJiUAlVi 


J 

J 


1 cnn n 
loUU.U 


10 ift_Ln iQ 

ly. Id it u. lo 


zUUy VJCt l.io 


A CiQA 1 

4.yo41 


T>"DOi\/nr>T' 
riXKjNir 1 - Z 


V 


ZvSOzU.U 


> zl.o4 


zUUy yJct L.l 1 


A OQO(^ 

4.yoyD 


Jr^ixUMJr^ i - 


TJ 

r> 


/1 /I n n 
z4z4U.U 


1 77 _L n 1 7 
zl. (I ± U.17 


zUUy VJCt i.zo 
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zUUy VJCt i.D4 


/I ni c; 
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U 
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zUUy UCt Z.U4 


QC7Q 
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zU. ( z it U. 11 
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Z04UU.U 


> ZZ. / 
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r> 
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1 


1 c^n-i n n 
10U4U.U 


on c;n 1 n 1 
zU.OU it U.lO 
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4z4U.U 
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zUUy UCt O.U4 


D.OOi ( 
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4UUU.U 
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Jr^ixWiVlJr^ ± -"y: 


E> 

-ft 
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2009 Oct 5.17 


9.0014 


PROMPT-5 


7 


13280.0 


> 21.24 


2009 Oct 5.49 


9.3449 


UVOT 


U 


4640.7 


22.07 ±0.16 


2009 Oct 6.38 


10.2257 


UVOT 


U 


4392.5 


22.84 ±0.29 


2009 Oct 7.17 


11.0096 


PROMPT-4 


R 


27840.0 


21.98 ±0.22 


2009 Oct 7.21 


11.0377 


PROMPT-5 


I 


23680.0 


> 22.04 


2009 Oct 7.42 


11.2777 


UVOT 


U 


5096.2 


22.22 ±0.17 


2009 Oct 8.32 


12.1742 


UVOT 


U 


4573.2 


22.94 ±0.33 


2009 Oct 19.05 


22.90 


Gemini-S / CMOS 


9' 


600.0 


24.56 ±0.11 


2009 Oct 19.05 


22.90 


Gemini-S / CMOS 


r' 


600.0 


23.97 ±0.11 


2009 Oct 19.05 


22.90 


Gemini-S / GMOS 


i' 


600.0 


23.67 ±0.12 



UT at beginning of exposure. 
^ Time from midpoint of exposure to FermvGBM trigger. 

^ Reported magnitudes have not been corrected for Galactic extinction {E{B — V) — 0.024 mag; ISchlegel et al.||199^ . Observations in the 3'-, 
r'-, and z'-bands are reported on the AB magnitude system ]Oke Gunn|fl983l ). All other filters are referenced to Vega. 



